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61 . THERMAL ANALYSIS OF CONNECTORS 
by 

S. Levy and L. H. Tomlinson 

61.0 Summary 

This Section has been divided into four parts: Section 61.1, Fluid 
Heat Transfer Coefficients; Section 61.2, Radial Heat Flow Through a Surface 
Film into a Circular Flange; Section 61.3, Temperature Transients in Bolts; 
and Section 61.4, Temperature Transient Analysis in Flange Assemblies. 

Section 61.1 is included because knowledge of the fluid heat- transfer 
coefficient is necessary in order to evaluate the magnitude and duration of 
temperature gradients within a flange connection during the transient state 
after flow has started. In several examples, coefficients are calculated for 
a hot combustion gas, liquid oxygen in straight forced convection, and liquid 
oxygen in two-phase flow. Also a coefficient is calculated for a blanket of 
stagnant gas to show how it could be used to insulate the inner wall of a 
flange against sharp temperature changes. It may be stated that calculations 
of heat- transfer coefficients in single-phase flow are fairly simple and 
straight-forward as illustrated in the first two examples. In two-phase flow, 
however, this is not the case, and in many situations the coefficients may 
only be estimated or guessed at by people with considerable experience and 
knowledge in the field. 

Section 61.2 presents a thermal analysis of transient heat flow into 
a flange which takes into account the curvature of the flange and is the 
type of analysis necessary for the larger radius ratios (outer radius divided 
by inner radius. It is assumed that the flange temperature rise (or fall) can 
be adequately described by a function such as T = Ar n 4- B where A and B are 
functions of time and n is a constant. A and B are determined by the conditions 
that (1) the total change in heat content of the flange equals the integrated 
heat flow through the boundary layer film; and (2) the total change in heat 
content of the flange outside the mid-radius circle equals the integrated heat 
flow through the mid-radius circle. The constant n is determined to satisfy 
the requirement of equal heat flow through the surface film and into the flange 
at the time of maximum temperature difference from inner to outer wall of the 
flange. It is assumed that the heat flows radially only. The film coefficient 
and the thermal properties of the flange material are considered constants. 
Analytical expressions are given for the maximum temperature difference and 
the time at which it occurs. The results are in a form which can readily 
be used in Section 48 to obtain distortions and stresses. Numerical results 
are presented for several examples. These show a greater temperature difference 
for the case of a stainless-steel flange than for an aluminum flange, as might 
be expected, The results are also compared with results using "exact" methods 
for a flat slab Of thickness equal to the flange thickness. Multiplying 
certain coefficients by a factor of 0.95 brought the agreement to within about 
a percent. 
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Section 61.3 presents some preliminary information on what the bolt 
temperature transients might be. Some temperature transient curves worked up 
for a particular flange are presented and are used also for the bolt which 
has a similar path length. Also included is a brief discussion of thermal 
resistance of metal- to-metal surface contacts along with curves summarizing 
some of the most applicable data. These data show that for low bolt pressures, 
and especially for the steels, the thermal conductivity may run from a few 
hundred to a few thousand Btu/hr.ft? °F. At high bolt loading and for the 
softer metal, aluminum, the thermal conductivity may reach many thousands of 
Btu/hr.ftr °F and the joint may for all practical purposes be considered to 
be bridged by conducting metal. 

Section 61.4 presents the analysis of temperature transients in flange 
assemblies by the use of a digital computer. Preliminary investigation shows 
that the low-profile flanges (with radius ratios close to 1) used commonly in 
the missile industry may be regarded as sections of flat plate for purposes 
of analysis. The permissible radius ratios are as high as 1.4. Analysis of 
a simple model system representing a flat plate on the digital computer showed 
very good agreement with previously known flat-plate data. Equations are 
presented for the analysis of much more complex model systems representing a 
specific flange assembly. The computer solution for this case is given. It 
shows good agreement with the appropriate results obtained in Section 61.3. 
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61.1 


Fluid Heat Transfer Coefficients (by L.H. Tomlinson) 


61.1.1 Introduction 

The fluid heat transfer coefficient at the inner wall of a flange 
^ill have an important effect on the magnitude and duration of thermal 
transients within the flange. Very high coefficients (in the thousands of 
Btu/hr. ft. °F) result in the inner surface driving toward the fluid temp- 
erature very quickly. In this case, thermal constants of the flange material 
exercise most of the control over the thermal history. Lower heat transfer 
coefficients which may sometimes be encountered will result in reduction of 
the temperature transients and gradients, and the magnitude of the coeffi- 
cient will be of considerable importance. It should be noted at this time 
that all work to date has considered the outside of the flange insulated. 

Heat flow from the pipe wall into the flange will be small until a combina- 
tion of very low fluid heat transfer coefficients and high flange conductive 
ity is encountered. 

The thermal characteristics of the system may be best described 
through the use of two dimensionless numbers: 

2 

the Fourier Number, Kr/pcs 
and the Biot Number, hs/K 

where, K = thermal conductivity of flange material 
p = density of flange material 
c = specific heat of flange material 
s = radial thickness of flange 
T = time 

h = heat transfer coefficient of the fluid at inner surface of the 
flange 

It will be shown elsewhere how these numbers are related to the 
thermal transients in the flanges. For the remainder of this section we will 
confine our discussion to methods of evaluating the h in the Biot Number 
above . 

The heat trhnsfer equation with the widest application to the fluids 
encountered in missile work is the Dittus-Boelter equation for flow in pipes. 
In various standard heat transfer texts (Refs. 1-3), this equation may be 
found with several slightly different coefficients. For calculations done 
here, the following forms of the Dittus-Boelter equation are used (Ref. 2, 
p. 394) : 

Nu = 0.023 (Re)°* 8 (Pr)°‘ ^ for heating the fluid (1) 

Nu = 0.023 (Re)°* 8 (Pr)°’ 3 for cooling the fluid (2) 


The above dimensionless numbers may be written as follows: 


Nu = hD/K f Nusselt's number 
Re = DG/n Reynolds' number 
Pr = cn/K^ Prandtl's number 


(3) 

(4) 

(5) 
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where h = heat transfer coefficient of the fluid 
D » diameter of pipe 

= thermal conductivity of fluid 
G * mass flow density 
1 1 » fluid viscosity 
c * heat capacity of fluid 

Other equations will be introduced where required 
now be calculated are listed as follows: 

1. Hot combustion gas flowing in 2 in. I.D. pipe 

2. Liquid oxygen flowing in 20 in. I.D, pipe 

3. Nucleate boiling of liquid oxygen 

4. Stagnant film of oxygen gas 
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Typical set of 
consistent units 

(Btu/hr.ft? °F) 
(ft.) 

(Btu/hr. ft. °F) 

(lb. /hr. ft?) 

(lb. /hr. ft.) 

Btu/lb. °F 

Examples that will 
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61,1,2 Sample Calculations of Heat Transfer Coefficients 


Example 1 . Heat transfer coefficient for hot combustion gas flowing 
in 2 in. I.D. pipe. 

Given data: 10000 lb. /hr. flow rate 

1700°F gas temp, 

940 psi pressure 

The following were obtained from Ref. 3, assuming that values for air 
are close enough at the low temperature of the gases: 

p - 3.1 x 10 5 lb. /sec. ft. or 0.111 lb. /hr. ft. 
c =0.28 Btu/lb. °F 

= 0.044 Btu/hr. ft. °F 


It may be calculated that: 

10000 x 144 

W — — — — — — — 


= 4.58 x 10 5 lb. /hr. ft? 


D = 12 " 0,167 ft ' 

Equation (2) will be used for this case 

Nu = 0.023(Re)° ,8 (Pr) 03 

0.167 x 458000 , „„ ,„5 

Re = DG/p = = 6.89 x 10 

„ 0.28 x 0 . Ill „ 

Pr = Cp/K f = 0^444 =0.705 

Nu = hD/K f = 0. 0 23 ( 689000) °‘ 8 (0 . 705)°" 3 


h = 


0.044 x 0.023(46600) (0.90) 
0. 167 

254 Btu/hr, ft? °F 


Example 2 . Heat transfer coefficient for liquid oxygen flowing in 
20 in. I.D. pipe at a flow rate of 3917 lb. /sec. 


The following were obtained from Ref. 5: 

p = 0.36 lb. /ft. hr. 
c = 0.4 Btu/lb. °F 

= 0.08 Btu/hr. ft. °F 
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It may be calculated that: 

3917 x 3600 x 144 
G 100* 


= 6.45 x 10 6 lb. /hr. ft? 


D = ■—= 1.67 ft. 

Equation (1) will be used for this case, 

Nu = O.OZSCReJ^CPr) 0 * 4 

1.67 x 6.45 x 10 6 „„ ,„6 

Re = DG/p '~0~36 = 30 x 10 

0.4 x 0.36 , „ 

Pr = c|j./K f = o708 1,8 

Nu = hD/K f = 0.023(30 x 10 6 )°' 8 ( 1. 8° A 
, 0.08 x 0.023 x 9.6 x IQ 5 x 1.265 


= 1340 Btu/hr . ft? °F 


Example 3 . Heat transfer coefficient for nucleate boiling of liquid 
oxygen flowing at 3917 lb. /sec. in a 20 in. I.D. pipe. Pressure is low at 
this point (100 psl) . 

First the burnout heat flux for pool boiling (no flow) is calculated by 
the following equation of Zuber (Ref. 3, p. 349). 


Q/A = 


„ - , ,0.25 0.5 

"V. / p e +p v| 


The terms are defined and evaluated as follows, the information coming 
from Ref. 5. 


p, = 72 lb. /ft. 


liquid density 


p =0.25 lb. /ft; vapor density 

v 

h = 92 Btu/lb. heat of vaporization 

fg 

-3 2 

a s 35.3 3 ; 10 lb. /sec. surface tension 


g = 32.2 ft. /sec. 


gravitational constant 


Substituting these values in Eq. 6 


n x 92 x 0.25 35.3 x 10 x 32,2 x 72 


neglecting final term 
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Q/A = 18 Btu/ft, sec. 


65,000 Btu/ft. hr. 


i 

burnout heat flux 


Now that the burnout flux is known, a heat transfer coefficient for pool 
boiling is calculated by use of the following equation of Labuntsov- (Ref. 4) 


where 


h_ 

Kf 


cP . 


aT 


(p H - ) ‘ 
v fg 


°- 125 <V 


0 . 


%r) : 


N. 


Pj(Q/A) caTg 

Re n(p H ) 3 J 
N v fg 


(7) 

( 8 ) 


The additional terms needed for the above equations are defined and 
evaluated as follows, the information coming from Ref. 5 


K 




c 

J 

T 

£ 

P 

a 


=0,08 Btu/hr . ft. F 
=0.4 Btu/lb. °F 
= 778 ft. lb. /min. Btu 
= 205°R 

= 0,36 lb. /ft. hr. 

= 1.1 x 10 -3 lb. /ft. 

0.4 x 0.36 


Pr = C|a/Kg = 


thermal conductance of liquid oxygen 
thermal capacity of liquid oxygen 
mechanical equivalent of heat 
saturation temperature at 100 psi 
viscosity of liquid 
surface tension 


0.08 


= 1.8 


Evaluating equation (8) we get 

72 2 (65,000)0.4 x 1.1 x 10~ 3 x 162 
0.36(0.25x92) 778 

= 7.05 

Equation (7) may now be evaluated. 

1 

h.0.4 x 72 x 1.1 x 10' 3 x 162 . „ 125( 7-05) 0. 65 ( 1>g) 3 

0.08(0.25x92) 778 

h = 3230 Btu/hr. ft? °F pool boiling heat transfer coefficient. 
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The temperature difference between wall and saturated liquid may now be 
calculated 


T -T = Q/Ah = 
w x 


65000 

3230 


20 °F 


The next step is to compute a burnout heat flux for the flowing liquid 
oxygen by the following equation of Kutaleladze (Ref. 6) 


Q/A 

P H, 
v fg 


0.085 


°8<p / e-p v ) 


0.25 


V 


p i“ p v 

erg 


0 . 25-, 0 . 5 r 0.5 

1 + 0.057 — 

\ p vf 


V 1 * 


H 


fg 


(9) 


Two additional terms need to be defined and evaluated 

V = 25.3 ft. /sec. 

T = 205°F 

i 

Equation (9) is now evaluated 


velocity of liquid flow 
temperature of liquid 


Q/A 

,025 x 92 


0.25 


^025 

|_35. 3 x 10' 3 x 32.2 x 72 


25.3 

1 72 \°- 25 l 

0.5 

1+0.057| 

( 72 \°* 5 0.4(0) 

-3 A 


0.25 92 


[ 35, 3 x 10 x 32. 2/ 



i 1 J 


= 0.085 


Q/A = 9,700 Btu/sec. ft? °F 

2 

= 350,000 Ftu/hr . ft? °F 


Now Fig 61.0 is constructed as shown with a line of slope equal to 3 
passing through the pool boiling flux and (T w ~T g ) . Then the new burngut flux 
just calculated is entered on the curve and a new ( T W " T S ) value of 39 F 
obtained. This graphical method was described to the author in a personal 
communication from N. Zuber who is now in the Advanced Technology Laboratories 
of General Electric Co, 


Now h = 


Q/A 


<VV 

350,000 

39 


= 9,000 Btu/hr . 


ft? °F 


heat transfer coefficient for 
nucleate boiling in the flowing 
oxygen stream. 
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Example 4 . Heat transfer coefficient for a stagnant film of oxygen 
vapor at low temperature. 

The main reason for making this simple little calcuation is to show how 
low the heat transfer coefficient might go and to suggest a possible use 
for this. 

Assume the vapor film to be at about -250F and the thickness o to be 
1/16 in. Then Ref. 5 gives a value for K f of 0.0061 Btu/hr. ft. F, 

Then: h = K /s 

= 0.0061 x 12 x 16 
= 1.15 Btu/hr. ft? °F 

The actual value will no doubt be several times higher due to natural convec- 
tion or flow leakage, but this is still far lower than any of the normally 
occurring fluid coefficients and could be applied to certain flanges in 
order to reduce the magnitude of the temperature transients in them. The 
method of application might be to place a thin metal sleeve liner in the 
pipe at each flange location so that vapor might be trapped behind it. This 
scheme would be considerably less effective for hot gases due to heat radia- 
tion and to the higher conductivity of a hot gas. 
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61.1.2 Discussion 


From the above examples and other calculations made in the past, one 
may devise a chart illustrating the approximate ranges into which the 
coefficients may Be expected to fall. 


Flow condition 


Gas flow 


Heat transfer^coefficient 
(Btu/hr ft °F) 

100-500 


Liquid flow 

Straight convection 

Film boiling 

Nucleate boiling 


200-2000 

100-300 

1000 - 10,000 


The above ranges of coefficients are higher than would be expected in 
normal practice and most examples will perhaps lie on the high end of these 
ranges. These high heat transfer coefficients are commonly encountered in 
the missile field because fluid transfer is being carried out very rapidly in 
minimum- we fght equipment with flow friction losses being of secondary 
importance . 

In actual practice it is possible that in the initial transient, all 
regimes of heat transfer coefficients may be passed through by cryogenic 
fluids at lower pressures. Initially, with high temperature differences, 
film boiling may exist. Then as temperature differences decrease some 
nucleate boiling may take place, only to be followed by straight liquid 
convection as inner wall temperatures fall to near the fluid temperature. At 
higher pressures, forced convection may be the only mode of heat transfer 
encountered. The boundaries for each of these modes will depend upon: (1) 

the fluid properties; (2) the pressure and flow conditions; (3) the material 
and temperature of the flange; and (4) the conditions just upstream of the 
flange . 

11 
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61.1,4 Conclusion 


In conclusion, two main statements should be made on the subject of 
calculation of heat transfer coefficients, 

1. Calculations are straightforward for the case of single phase flow 
and use of the Dittus-Boelter equation or its equivalent will be satis- 
factory for practically all the fluids encountered (not liquid metals) . 

2. Calculations are not so simple for two-phase flow, and in some cases 
exact coefficients cannot be obtained without experimental data. This field 
is still being studied and data obtained and correlated. It would be wise 

to consult an authority on two-phase flow or at least have on hand authori- 
tative information covering the situation in question. 
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61.2.2 Nomenclature 


T 

A, B,n 

r 

a 

b 

t 

c 7 
h 


temperature rise (°F) . (Minus sign indicates a fall) 
terms in T = A(r/a) n + B 
radius (ft) 

inner radius of flange (ft) 
outer radius of flange (ft) 
time (hr) 

heat content per unit volume of flange material per 
F temperature rise (BTU/cu ft/°F) 

film coefficient (BTU/hr/sq ft/°F) 


k thermal conductivity (BTU/hr/sq ft/(°F/ft)) 

T inner wall temperature rise (°F) 

a ^ 

outer wall temperature rise ( F) 


T 

^i> C 2 9 C 5 


fluid temperature difference from initial flange 
temperature ( F) 

* m l t m 2 t 

terms in B/T = C^e + 


CL and C, 
3 4 

m l 9 m 2 


^ m_ t m t 

terms in A/T = CLe + C,e 
3 4 


constants in above equations 
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61.2.3 


Results 


It is shown in Section 61.2.6 that the temperature in the flange is 
given by: 

r 

T/T = + C 3 (r/a) exp(m t) + ; + C^(r/a) expCu^t) 


c M m 2 m l m 2 f 1 Wcy| L b 

1 \ m -m, i n+2 nk 

l 1 / \ 21/t M 


(n+2') _ ^ a+b^ 

W 


(n+2) 


( 0 . 95 ) 


C 9 = -C.-l (see Fig. 61.3) 


c J(!«-|/c2l| r ~ 1^1 

3 |K‘ m l | 2nk f»+b\ n 


C 4 = -C 3 (see Fig ' 61 ,3 ) 


(0.95) 


n = log 


1 4a+4b 


|^/log (see Fig. 61.1) 


and and m 2 (see Fig. 61.2) satisfy the equation 


C c m + C.m + C_ = 0 
-> o / 


[ b (n+2) . a <„+2>] (b 2. a 2 

(n+2) " (n+2) 


b 2 -a 2 ) L b 


(n+2) /a+b 


m 


(n+2) 


2 a+b 


ah n 
— a - 

cy 


(■,2 2 "I / a+b 1 n 

l b - a J — 


2ah 

cy(n+2) 


(n+2) a+b 


(n+2) 


[> 2 - ml 


2ah| kn 


cy Icy 


2 a+b 
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The maximum value of A (see Fig 0 61,6) is valuable in computing stress. It 
is given by 



18 
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,1.2.4 Examples 


Example 1 : Consider a flange joint in a pipe (stainless steel) which is 

suddenly exposed to the flow of liquid oxygen at a temperature 300 F less than 
the initial temperature of the pipe. ThenT* = -300 F. The inner flange radius 
a = 3 in. s 0.250 ft. The outer flange radius b = 4 in. = 0.3333 ft. £et the 
value of cy = 52.4 BTU/cu. ft./°F, the value of k = 12.4 BTU/hr./f t . / ( F/ft.) 
and the value of h = 500 BTU/hr./sq. ft./°F. What is the maximum temperature 
difference in the flange? 


From Eq. 1 or Fig. 61.1, n = -4.1 


b. 

c. 

d. 


From Eq. 2 or Fig. 61.2, m = -604, m 2 = -42.9 

From Eq. 4 or Fig. 61.4, t. N -0. 00471 hrs . = 1.70 sec. 

(max, diff.) 


From Eq. 5 or Fig. 61.5, (T -T.) 

’ a b 


- -171 F 


max 


e. From Eq. 1 or Fig. 61.3, ^ = 0.519 C 2 = -1.519 

C 0 =*1.015 C. = 1.015 
3 4 

f. From Eq. 1 at the time given by (c) and using Eq. 4, we get 


(T ) = - 172 F (T ) 

a (max, diff.) (max, diff.) 


= -1.2 F 


g- 

h. 


From Eq. 6 or Fig. 6 1.6 (A) = -250 F 

^ max 


It is of interest that at the time of maximum temperature differ- 
ence, exp (m^t) = 0.0579 and exp(m 2 t) * 0.8169. 


Example 2 : In this example we consider a larger flange with material 

having greater conductivity (aluminum). Inner radius a = 10 in. = 0.8333 ft., 
outer radius b = 11.3125 in. = 0.9427 ft., cy = 36.2 BTU/cu. ft./ F, 
k = 118 BTU/hr./sq. f t . / (°F/f t . ) , and h = 500 BTU/hr./sq. ft./ F. What 
is the maximum ratio of the temperature difference to the fluid temperature? 


a. Proceeding as was done in Example 1, we find n = -10.5, giving 
m 2 = -102.0 and m^ = -2425. 

b. Substituting these values in Eq. 4 or using Fig. 61.4, givefe 

t, ,. = 4.9 secs. (.00136 hrs.) 

(max. diff.) 

c. Substituting in Eq. 5 or using Fig. 61.5, gives (T -T ) = 0.176T* 

Si D 

max 


JO 
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61,2,5 Discussion 


An approximate method of calculating thermal transients in pipe 
flanges has been presented. It takes account of the flange curvature and 
in this respect adds additional information to "exact 11 methods that consider 
the flange to be flat. A correction factor of about 5 percent has been 
applied to some coefficients to bring the results presented within about 1 
percent of "fitting in" with the flat-plate results. Curvature is shown to 
be a small factor. No account has been taken of axial flow of heat from the 
hub into the flange. In general, it is felt that such a heat flow would in- 
crease the effective value of h and thus result in somewhat greater temp- 
erature differences than would otherwise be reached. 

The examples indicate that substantial temperature differences develop 
in flanges having low thermal conductivity in the presence of a high film 
coefficient. Flange thickness increases also increase the temperature dif- 
ference. For the examples presented, the time required to achieve the 
maximum temperature difference was only a few seconds. This indicates that 
thermal transients in pipe flanges occur quickly. 

In the design of pipe flanges, it might be possible to modify con- 
siderably the transient behavior by including an insulating sleeve at the 
flange or by modifications in the flow pattern at the flange. Such schemes 
of thermal design warrant further study. 


20 
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61.2.6 Appendix: Thermal Analysis 

We assume that the flange temperature rise can be adequately described 
by 

T = A(r/a) n + B (14 

where A and B are functions of time and n will be subsequently selected for 
fit. We can determine A and B by the conditions: 

a. The total rise in heat content of the flange equals the net heat 
flow through the boundary layer; and 

b. The total rise in heat content of the outer half of the flange 
(outside of mid-radius) equals the net heat flow through the 
mid-radius layer. 


These conditions give 


ah |T*-T 




rTcydr 


rTc7dr 


— k ^ 

2 I dr lr= a+b 

' ' 1 2 


Substituting Eq. (1A) into (2A) and (3A) and integrating 


$ (T*-A-B) - 


/ 2 2 1 

1 dB 

lb -a 

l 2 1 

dt + [ 


(n+2) (n+2) 

f " cL 

(n+2)a n 


_ dA 
n \ dt 


b (n+2) (n+2) 


log — as n 
& a 


nkj |a+b 
ic 7 2a 


n b 
A = - 


2 |a±b| 2 l , (n+2) f a+b | (n+2) - 

"1 2 I dB "1 2 I dA 

2 dt / n J dt 

J t (n+2) a J 


Note : 


(n+2) | a+b' 


(n+2) 


r 2b 

iog as n 


61-19 



Differentiating Eqs. (4A) and (5A) with respect to time, and using £he 2 
resulting equations together with Eq. (5A) to eliminate dB/dt and d B/dt 
gives 

Sft +C 6 f tC 7=" (6A) 

dt 

where C c , C,, and C_ are given with Eq. (2). 

00 / 

Equations such as (6A) have solutions in the form 

A/T* = C 3 exp (m^t) + exp(m 2 t) (7A) 

where and are arbitrary constants and m^ and m 2 are roots of 

C c m 2 + C,m + C-, = 0 ( 8A ) 

J o / 


Equation (2) is the same as Eq. (8A) . 


Substituting the value of A given by Eq. (7A) into (5A) and integrating 
gives 


*- c io 


-C 3 exp 


/ 

nk 

a+bl n 

2 

1 

f_2_] 

' b (n+2). ( a±b ) (l ' +2) 

\ 

C 7^ 

[2a j 

l 2 -m 2 \ 

T 

| n+2j 

a"b 2 -(^) 2 a" 

I 


f 


-C, exp (m^t)v 


nk 

c/m n 




(9A) 
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We know that initially T is zero for all values of r. Substituting (7A) and 
(9A) into (1A) and setting T = 0 and t = 0 gives 


0 = - C 


nk 


10 3/|c7m^ 


C 4 


nk 


cm. 


a+b 



2 


2a 

i 

P 

I 

a* 

ro 

1 

a+b 

2 

2 

/ a+bj 

n 

r 

2 

1 

2a 


b 2 - 

a+b 

j"2 




2 




(10A) 


2 

n+2 j 


b (n+2) .(i+b) 


(n+2) 


l\ 


A 2.(a+b)V 


o = c„ + c. 

3 4 


(11A) 


We know also that when t is very large, T = T . We will find that nu and m2 
are negative, so that when t is large, exp(mt) is zero. With (1A) then 

c 10 = X < 12A > 

Combining Eqs. (10A), (11A) and (12A) 


0 = 1 - C 


a+b 
3| 2a 


n 


L» 2 -(^)' 


mk 

C7j 


1 m l m 2, 


(13A) 


giving 




CL =i 
3 


c 7 
2nk 


b 2 -^) 

m n 


= -C, 


(14A) 


Substituting (7A), (9A), (12A) and (14A) into (1A) gives Eq. (1), except 
for the factor (0.95). This factor was applied to force the solution to fit 
in somewhat better with Ref. 9.' Without this factor, the curves for 
b/a = 1.1 and 1.3 in Fig. 61. 5 would have fallen above the curve for 
b/a = 1.0, taken from Ref. 9. Applying a f Actor of 0.95 to and 

caused the curves for b/a = 1.1 and 1.3 to be lowered as shown in Fig. 61. 5 
In this position they "fall in" with the flat curve b/a = 1. Since 
; coefficients C 1 and C also deal with the transient heat flow phenomena, 
it was felt that the 0.95 factor should be applied to them also. Apparently 
requiring the temperature distribution to arbitrarily fit the form of Eq. (1A) 
with only two parameters, A and B, has introduced this 5 percent error. The 
solution has shown that changing b/a from 1.0 to 1.3 has only a small effect 
on the answers. This raises the question of using the results in Ref. 9 
directly, neglecting curvature effects. This was not done, since a solution 
in the form of Eq. (1A) is particularly useful in subsequent stress and dis- 
tortion computations, Section 42, 
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The temperature difference between the inner and outer flange surfaces 
is obtained by evaluating T and T, at r = a and b, respectively, from 
Eq. (1) and forming their difference. The result is given in Eq. (3). To 
obtain the time at which the difference is a maximum, the derivative is set 
to zero. The resulting solution for the time of the maximum temperature 
difference is given in Eq. (4). Substituting the resuUinto Eg. (3) gives 
Eq. (5). From Eq. (1A) it is evident that A = (T a ~T^)[ a / ( a "b )]• Using this 
result in Eq. (5) gives Eq. (6) 


We will now select n to give equality of heat flow through the film and 
at the inner face of the flange at the time given by Eq. (4). In equation 
form this is expressed by 


h(T*-T ) = -k| 
a 


dT 

dr 


when t * t 


r=a 


(max, diff^ 


(15A) 


Substituting from Eq. (1) into this equation (neglecting the 0.95 factor for 
this purpose) gives 


m 0 C 0 + — C = nu C 4- — C Q 
2 8 m^ 9 18 m 2 9 


where 


C 8 = 


. nk n 

1 a - 

ah 


2 

n+2 


L (n+2) fa+b\ 
, b ~1~*2~) 




(16A) 


[* 2 ] 



Equation (16A) can be put in the form 


n^'d* m^ 

tT, 



Likewise from Eq. 2 we know that 


(17A) 


I m 2 + m^ 
m l m 2 ] 



(18A) 


Combining Eq. (17A) and (18A) and solving for n gives 


log 

[a + 3b ] 

14a + 4b 

log 

a + b| 

hc”J 


This is the value of n given with Eq. (1) and in Fig. 61.1 









k 


Fig. 61.4 Dimensionless time of maximum temperature difference 
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61 . 3 Temperature Transients in Bolts (by L,H. Tomlinson.) 

61*3.1 Introdu c tion 

Knowledge of the temperature transients in the connector bolts is 
necessary in order to determine the changes in stress and strain in the bolt 
and the overall effect on the connector. In the case of a sharp and large 
heating transient in the connector, it may be possible to actually over- 
stress' the bolt and put a permanent set in it before it has had a chance to 
come up to temperature. 

In the case of a sharp and large cooling transient, it may be possible 
to release all bolt loading on the connector before the bolt gets a chance to 
come to temperature equilibrium. 

Determination of the actual temperature transient in the bolt in relation 
to the fluid transient is quite complex and requires considerable mathe- 
matical effort, preferably on a computer. Efforts of this nature will not be 
covered here at this time, however. Instead, some simplified methods of 
looking at the problem will be discussed. 

Considering the bolt shown in Fig. 61. 7 one assumption that may be made 
is that there is no heat transfer between the bolt and the connector other 
than by contact with the bolt head and nut. At room temperatures and below, 
this is very nearly true^ as the convection coefficient around the bolt will 
be only 2 or 3 Btu/hr.ftT °F and the radiation coefficient will be negligible. 
At high temperatures the coefficients can be considerably higher and might 
run as high as 20-40 Btu/hr.ft. F for connectors carrying hot combustion 
gases . 

Accepting the simplifying assumption of only lengthwise heat flow in the 
bolt, one may look upon the bolt as merely a small section of an infinite 
flat plate and apply some of the already available thermal transient infor- 
mation to certain cases. In order to get a rough picture of what the thermal 
transients may look like, Fig.61.8 was included. This graph was drawn 
originally to show the temperature transients occurring in flange SK20-1286A 
(sketched in Fig. 61. 7) with liquid oxygen as the fluid and a heat transfer 
coefficient at the fluid-wall interface of 2000 Btu/hr.ft. °F. This curve 
was drawn from flat plate data in Ref. 1 and Ref. 2 using dimensionless 
numbers containing the thermal constants of stainless steel and aluminuim 
the dimensions of Fig. 61. 7 and the above mentioned h of 2000 Btu/hr.ft. °F. 
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CROSS SECTION OF 
BOLT AND FLANGE CONFIGURATION 
SIMILAR TO SK 20-1286A 



Scale 1 to 1 























13. 2 Thermal Resistance at Contact Points 


Before we can even use the curves of Fig. 61.8 as approximations for some 
cases, we must discuss the matter of the contact resistance at the bolt head 
and nut. Because contact resistance is quite dependent on surface finish as 
well as on material and the pressure on the joint, only general approxi- 
mations can be given here. Fig. 61.9 shows some curves (from Refs. 10 an ) 
that should apply reasonably well. These data apply to surface roughnesses 
in the range of 50-100 microinches, which is equivalant to a fairly ordinary 
machining job that might be done on a flange or bolt. Certainly it would be 
unreasonable to expect contact surfaces on bolts and flanges to possess a 
high degree of flatness and a smoothness on the order of 10 microinches or 
less . 


It will be noted that the hardest pair of surfaces, steel on steel, 
exhibit the lowest thermal conductance at the contact area. Both stainless 
steel on aluminum and aluminum on aluminum show much higher thermal 
conductances due to the much greater tendency of the softer materials to 
deform locally at their scattered points of contact and increase the area of 
actual physical contact as pressure is increased. Sometimes as the pressure 
raised very high, very considerable amounts of plastic flow occur, even 
in steel, and the resulting high thermal conduction is retained down to very 
low pressures. 

A point to note, especially at lower conductances, is that at the nut 
end of the bolt there are two resistances in series. 
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>1.3.3 Examples 


It will be shown how flat-plate data may be used for getting rough esti- 
mates of temperature transients in symetrically heated bolts. It should be 
kept in mind while using flat-plate data from Ref. 1, pp, 33-34^ and Ref. 9 or while 
referring to Fig. 61„8 that the following terms are equivalents: 

Flange outer surface=Bolt center=Flat plate midplane 

Flange inner surface=Bolt ends=Flat plate outer surfaces 

Flange midplane=l/4 way thru bolt=l/4 way thru flat plate 

The reasons for these terms being equivalent may be seen by considering 
the following facts: 

1. Heat flow into the infinite flat plate is from both sides and at 
the midplane there is no heat flow. 

2. Heat flow into the flange is from one side but the outside is 
considered insulated and there is essentially no heat flow thru the outer 
boundary. 

3. Heat flow into the bolts is limited arbitrarily in our cases to 
symetrical flow from both ends. Thus there is no heat flow thru the center 
of the bolt. 

It can be seen therefore that the heat flow through the flange is 
equivalent to heat flow through one-half the thickness of the flat plate or 
through one half the length of the bolt. 

(a) First case - Stainless steel bolt as per Fig. 61.7 

Stainless steel flanges as per Fig. 61.7 

Assume a contact surface pressure of about 7,000 psi, producing thermal 
surface conductance of about 2000 Btu/hr. ft. °F. For a first look at the 
bolt, we could assume that the temperature of the flange had changed 
instantly. Than both halves of the bolt would exhibit temperature transients 
as in Fig. 61,8 with the center of the bolt being equivalent to the curve 
labeled "outer surface"- It is evident that at this rate, 100 seconds will 
elapse before the center of the bolt comes down even halfway in temperature 
and much longer before the bolt as a whole comes close to the final temper- 
ature . 

Since the assumption of instantaneous temperature change in the flange 
is certainly not correct, it is well to take a look at what actually happens. 

It will be noted from Fig. 61.8 that the point at which the bolt contacts the 
flange (about 3/4 of the way out) really follows a temperature curve rather 
closer to the outer surface than the midplane. This, then, is the initial 
temperature applied, thru the contact resistance, to the bolt head or nut. 

Obviously the temperatures in the bolt will take even longer than previously 
estimated to react to the initial fluid temperature change, A good estimate 
might be 2.5 minutes for the center of the bolt to come down halfway in 
temperature . 
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(b) Second case - Aluminum bolt as per Fig. 61.7 

- Aluminum flange as per Fig. 61.7 

Assume a modest contact surfac«> pressure to produce a thermal surface 
conductance of about 2000 Btu/hr.ft. ®F. • Actually, though, higher thermal 
conductivities are not going to alter greatly the answers arrived at here. 
Again making a temporary assumption that the flange temperature changes 
instantly, it becomes evident from Fig.61.8 that the center of the bolt will 
be halfway down in temperature in about 14 seconds. If the actual rate of 
change in flange temperature is taken into account, it is estimated that it 
will take about 25 seconds for the center of the bolt to drop halfway in 
temperature . 

Other more complex cases, where flanges of different materials are used, 
will not be discussed here, as even rough estimates of temperatures in the 
bolts will be difficult and hardly worthwhile. 
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61. 4 Temperature Transient Analysis in Flange Assemblies (by L.H. Tomlinson) 

61,4,1 Preliminary Computer Efforts 

Sections 61,2 and 61,3 present information on temperature transients 
for plain flanges and for flange bolts in certain simplified cases. There 
is a need for some knowledge of the thermal transients in more complex assem- 
blies, however. This could cover combinations of flanges, bolts, spacers, 
and sometimes, gaskets of various materials in the same assembly. 

One question that needed to be answered was whether or not the flanges 
could be considered as sections of flat plates due to their relatively low 
profiles or low radius ratios. Ref, 12 gives data on thermal transients in 
cylinders, and some of it lies in the range in which we are interested. Fig. 
61,10 presents a curve plotted from flat plate data given in Refs. 1, pp. 
33-34, and Ref. 9, Data from Ref, 12 for radius ratios of 1.3 and 1,4 are 
plotted on Fig, 61,10 and it will be noted that the points fall very close to 
the curve for flat plate transients. The conclusion to be drawn is that we 
may consider our low profile flanges of radius ratio 1.1 to 1.3 to be sections 
of flat plate. 

As a test, a section of flange or plate was divided into nodes as 
shown in Fig. 61.11 and the dimensionless heat transfer equations of Tables 
I and II were derived. A very brief description of their writing follows. A 
heat balance equation is written for the node "a" as the generalized example 
following shows: 

A (t v -t ) + B (t -t ) + ^ = 0 

v b a' v c a' AT 

This equation is written with so-called "forward time steps" wherein all the 
temperatures t are known for the beginning of the time step and only one 
temperature t T of the node itself is unknown and therefore readily solved for. 
The other way of writing the equation with "backward time steps" is as 
follows: r 

a + b (t; t;> + c < £ a > = o 

Here there is only one known t and the rest unknown. This results in a 
matrix of unknowns to solve for each time step. This latter method is in- 
herently a much more accurate one but it is very time consuming to solve a 
matrix each time. The use of "forward time steps" avoids this difficulty 
but introduces possible error because of the assumption that 



AT 


is constant throughout the time interval. Since there is usually more or 
less curvature to the function the only way to keep error down is to make 
the time interval (AT) quite small. It was decided to try "forward time steps" 
since there could be considerable savings in computer time and in programming 
if the time steps did not have to be extremely small. The equations of Table 
I were made dimensionless by appropriate manipulation of the constants in- 
volved, Trial solutions were run on these equations using time steps of one- 
tenth of a second and printouts every second but these solutions proved 
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to be somewhat in error. Selection of a time step of one hundredth of 
a second with a printout of data every second gave results of quite satisfactory 
accuracy with only a small increase in computer time (to a little over 5 

minutes on the GE 225), 


-2 
. 11 

Division of thickness into nodes for step-by-step 
numerical solution 

It should be noted here that Ref. 13 contains a very good description of the 
methods being used here. Pages 19 and 20 of this reference have some especially 
pertinent comments concerning time steps. 
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TABLE I 


Equations for Simple Flat Plate Representation of Flange 
(written in completely dimensionless form) 


5. 


6 . 


7. 


R 1 = 


i + 


1 + A 


AAT 


1 + A 


R„ 


C(l+A ) 


AAT 


C (1+A ) 


+ Jgt-i 


R 4- bat _ 
R 2 + C R 3 


3. R_ = 


4. 


BAT 

E 


R 2 " 


DAT 
G R 3 


-[ 


RAT , J^T , 
~T + E ' 1 


BAT FAT . 

G G 


n .DAT 
R 3 + E R 4 


R 

R 4 + G 


if r 4 . [ m + »f . i] r 5+ ^i 
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TABLE II 


R 

n 

S 

R 

n 

AT = 
A = 
A = 

B = 
C = 

D = 
E - 

F = 
G = 

H = 
I = 

J = 
L = 

M = 

K » 
P = 
c = 
s = 
h = 


Definition of Terms Used in Equations of TABLE I 


t _t 

Known dimensionless temperature at beginning of time step(- — — — ) 

o f 
t - t ’ 

Unknown dimensionless temperature at end of time step (— — —^) 

o o 

time step in seconds 
40 


40 


hs 


10 


0.05 


pcs 


5.71 


0.15 


PCS 


0.2 


pcs 


0 2 £- 2 - 
u,z k 


0.2 


Pcs 


0.2 


Pcs 


thermal conductivity of flange material BTU/hr.ft. F 

3 


density of flange 

specific heat of flange material 

radial thickness of flange 


lb/ft' 
BTU/lb°F 
ft 


2o t 


heat transfer coefficient of fluid at BTU/hr.ft F 
surface 
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FIGURE 61. 12 


Bolt & Flange Assembly 
(cross section) 
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61,4.2 General Equations for Computer Solution of Flange Assembly. .Temp erature Tra nsients 


Fig. 61. 12 shows a combination of flanges, bolt and spacer each of which 
may be a different material. The nodal structure chosen is also shown. 

Heat balance equations were written for each node and the general equations 
of Tables III and IV obtained. These equations are non-dimensional only 
in temperature in order to make substitution of different materials and 
dimensions of flanges easier. 

It would be appropriate at this time to comment on why the flanges are 
left disconnected from the pipe in this analysis. Actually for most cases 
of cryogenic fluid flow the fluid heat transfer coefficient is so high that 
the pipe wall makes a rather ineffective fin for the flange. For those 
cases where the fluid heat transfer coefficient is not very high (perhaps 
under 500 BTU/hr.ft °F) an additional node will have to be introduced to 
represent a piece of pipe wall. 
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TABLE III 


Equations Written for Flange & Bolt Assembly of Figure 61.12 


z n AT "l „ r (z £i +z 2i> AI .] r ,hi* . 

— R f • — 'J 1 Sj I R 2 


Z 12 AT | r < Z 12 +Z 32 )AT , I . ,[ Z 32 AT . 

ST 1 ' % 2 S 2 3 


Z 23^ "^ Z 23 +Z 43 +Z 63^ A *' Z 43 A>r Z 63 AT 

-|H R 2 - [— X - l J R 3 \ + |77-J R 6 


Z 34 Af L f (Z 34 +Z 54 +Z 74> AT " , Z 74 AT . F Z 14-4 AI I, 

~J 3 ' L % l J *♦ L «4 J 5 L s 4 J 7 + — S 7 -J 


"z, ,-AT f(Z +Z )AT 1 "Zoc-AT 

-|5_ R . _4L_85 1 R + R 

s 5 4 S 5 5 L S 5 J 8 


Z 36 AT 1 „ R Z 36 +Z 76 +Z 11-5 )AT f] p76 AT Z 11-6 AT 

~ 3 ' L h _ 6 ptJ 7 nr J 11 


(Z 47 +Z 67 4 ' Z 87 +Z 10-7^ AT 


- 1 R 


Z, ..AT Z. _AT 

47 „ . _67 

S 7 R 4 + s7J6 


Z «7 AI Z 10-7 AI 

+ |-Ir Rr + “57^ R 10 

Z 58 AT 1 F Z 78 AT 1 r< Z 58 +Z 78 +Z 98> AI ,1 „ , „ 

^rJ R5+ L~J 7 'L h J 8 hrJ 9 

Z 89 AT 1 r P^lo-S^ .1. r Zl °-9 AT " r 

■vr® 1 ~ ' J 9 L s 9 J 10 

Z 7- 10 AT 1 „ JW*L F l - Z 7-10' t ' Z 9-10 t Z ll- 10- >AT ,1, 

~J 7 pFTj 9 1 s w ' J 10 

Z 11-I0 AI . Z 19- 10 AT ] 

s 10 11 l S Io J 19 
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TABLE IV 


Definition of Terms in TABLE III 

Z = U A 
n . n n 


fl 


'fl 




— H- 

h f K X 


, = - 6-11 
'&-11 1 /& ill 

h c *6 K n 


K 12 A 12 


'12 


7 


12 


' 9-10 


K 9- 10 A 9- 10 


X 


9-10 


K 23 A 23 


'23 


/ 


23 


' 10-11 


K 10 - 11 A 10-11 


7 


10-11 


K 34 A 34 


34 ^34 


' 11-12 


K ll-12 A ll- 12 


11-12 


"43 


K 45 A 45 


45 


K 12-13 A 12- 13 


'12-13 


77. 


13 


'36 


a 36 

1 ,^36 
h 6 K 36 


'f 13 


f 13 

ITT 


13 


h f K 13 


'47 


'58 


47 

i + hi 

h c K 47 
A 58 

1 758 

\ K 58 


'4-14 


'14-15* 


4-14 

nr 


h + K 
c 4 

K 14-15 A 14- 15 




15 


89 


'89 


rTTX 


hc + 


*9 


K„ 


K 15-16 A 15-16 


'15-16* 


77 


16 


"7-10 


h c K 7 K 10 


K 16- 17 A 16-17 


'16-17 


77 


17 
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TABLE IV 
(continued) 


Z 17- 18 


'18- 19 


'18-19 


K 17- 18 A 17- 18 


hi- 


18 


K 18- 19 A 18-19 


7 18- 


19 


lO- 19 

i iio-iT 

h c K 10- 19 


Further definitions 


n 


M 

n 

C 

n 

A 

nm 


h 


h 

i 


nm 


M C 
n n 


mass of node n 
specific heat of node n 

2 

cross section of node interface between n & m (ft ) 
fluid heat transfer coefficient 
contact heat transfer coefficient 
length of heat path in node n 
length of heat path in nodes n & m 


(BTU/hr ft 2 °F) 
(BTU/hr ft 2 °F) 
(ft) 

(ft) 


K = thermal conductivity in node n 


(BTU/hr f t°F ) 


n 
R 

n 

t 

C o = 

I 

R 

n 

AT = 


— — = temperature ratio at beginning of time stop 

t f t o 

nodal temperature at time T 
fluid temperature 
nodal temperature at time 0 
temperature ratio at end of time step 
length of time step 
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61,4.3 Computer Solution for Specific Flange Assembly 


The general equations previously presented have been solved for a 
specific case. This is for the test flange assembly shown in Figure 61.7 
where one flange is aluminum, the other is stainless steel, the spacer is 
stainless steel and the bolt is aluminum. Thus nodes 1 through 8 are stain- 
less and' nodes 9 through 19 are aluminum. The following physical char- 
acteristics are assumed for these materials: 

Aluminum 


K = 80 Bt-u/hr, ft.°F 

C = 0.211 Btu/lb, °F 

p = 0. 101 lb. /in? 

Stainless steel 

K = 10 Btu/hr. ft. °F 

C = 0.117 Btu/lb. °F 
p = 0.283 lb. /in? 


When the above constants and the physical dimensions are substi- 
tuted in the general equations of Tables III & IV, the equations of Table 
V are arrived at. 

The equations of Table V were solved by a digital computer using 
time steps of 0.01 sec. The program was run for 50 sec. with a printout 
of results every second. Figure 61.13 shows some of the most significant 
of these results plotted in terms of dimensionless temperature vs, time. 

The assumption is, of course, that the flange assembly is at temp- 
erature t Q when suddenly at time zero a fluid of temperature tf is passed 
over the inner surface. Temperatures then start towards tf and are given 
as a fraction of the total initial temperature difference t - t^. The 
bolt temperature B is just an arithmetic average of all the temperatures 
of the bolt nodes but it should be quite close to the true average temp- 
erature. It is of interest to compare some of the curves of this figure 
with those of Figure 61.8 which was prepared for just plain disconnected 
flanges of the same dimensions and materials. It is apparent that they 
show fair agreement. 
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Table V 


Equations for Test Flange Assembly of SK20-1286A 

1. 1.17 R1 =0.000853 R f + 1.167937 R 1 + 0.00121 R 2 

2. 1.17 R2 = 0.00121 R l + 1.16758 R £ + 0.00121 R 3 

3. 1.17 R g = 0.00121 R 2 + 1.167408 R 3 + 0.00121 R 4 + 0.000172 R g 

4. 1.17 R.' = 0.00121 R. + 1.166947 R, + 0.00121 R + 0.000172 R 

4 J 4 3 +0. 000461 R ' 

14 

5. 1.17 R^ = 0.00121 R 4 + 1.168618 R 5 + 0.000172 R g 

6. 0.219 R g = 0.000172 R 5 + 0.218048 R g + 0.000227 R ? + 0.000553 R n 

7. 0.219 R' = 0.000172 R. + 0.000227 R, + 0.217821 R + 0.000227 R„ 

7 4 6 + 0.060553 R in 8 

8. 0.219 Rl = 0.000172 R c + 0.000227 R^ + 0.218048 R„ + 0.000553 R 

8 5 / o y 

9. 0.753 R^ = 0.000553 R g + 0.742787 R g + 0.00966 R 1() 

10. 0.753 R^ 0 = 0.000553 R ? + 0.00966 R g + 0.732666 R^ + 0.00966 R u 

+ 0.000461 R 19 

11. 0.753 RJ = 0.000553 R g + 0.00966 R^ 4- 0.733127 R n + 0.00966 R 12 

12. 0.753 Rj 2 = 0.00966 R n + 0.73368 R 12 + 0.00966 R lg 

13. 0.753 R^ 3 = 0.00123 R f + 0.00966 R 12 + 0.74211 R lg 

14. 0.117 R' 4 = 0.000461 R 4 + 0.11620 + 0.000339 R 15 

15. 0.103 R’ = 0.000339 + 0.102429 R 1$ + 0.000232 R^ 

16. 0.103 R’ = 0.000232 R,. + 0.102536 R , + 0.000232 R 

16 15 Id 1/ 

17. 0.103 R’ = 0.000232 R., + 0.102536 R + 0.000232 R 

17 16 1/ to 

18. 0.103 R’ = 0.000232 R + 0.102429 R lg + 0.000339 R J9 

19. 0.141 R’ = 0.000461 R + 0.000339 R 1C + 0.14020 R 

19 1U to Ly 
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62. EFFECT OF PRESSURE SURGES 


by 

B. T. Fang 


62.0 Summary 

The pressure surges ("water hammer") in a pipeline system have been 
studied extensively by hydraulic engineers (Refs. 1, 2 and 3). Attention was 
mainly given to the severe case of sudden stop of liquid flow in a pipeline. 
For launch vehicles we are faced with the pressure surges in connection with 
the sudden start of flow. In this section the magnitudes and frequencies of 
these pressure surges are studied. The essential effects they have on the 
sealing of fluid connectors are 

(1) The higher pressure associated with these surges aggravates the 
already existing adverse effect of internal fluid pressure, 

(2) These pressure surges are sources of shock loading which may 
excite undesirable vibrations. 

(3) The results indicate that a proof pressure of twice the operating 
pressure would exceed any expected surge pressure. 
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62. 1 Water Hammer Effect as a Phenomenon of Interchange of Kinetic and 

Pressure Energy 

The water hammer effect is a phenomenon of the interchange of the 
kinetic and elastic (pressure) energy of the fluid. When the motion of the 
fluid is suddenly stopped, most of its kinetic energy is transformed into 
its elastic energy with a resulting increase in its pressure. Conversely, 
when the pressure of the fluid is suddenly released, most of its elastic 
energy is transformed into kinetic energy with a resulting increase in its 
velocity. It can be shown from an energy balance that the relation between 
the pressure change and velocity change is given by the simple formula 
(Ref. 1) 


Pressure Increment. Ap 
Velocity Reduction, -Av 



( 1 ) 


where p is the density and K the bulk modulus of the fluid. 

Equation (1) shows that the denser and the more incompressible the 
fluid, the greater is the magnitude of pressure surges. 
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62- 2 Transmission and Reflection of Pressure Signals 


Pressure signals are transmitted in a fluid with finite velocity - the 
velocity of sound. In a rigid pipe the velocity of sound of the fluid is 
(Ref. 1) 



For not-so-rigid pipes the velocity of sound is smaller. It can be calcu- 
lated from Eq. (2) with the bulk modulus K replaced by the apparent bulk 
modulus (Ref. 1), 

K' - K/(l + ~ K) (3) 

where r is the radius, t the thickness and E the Young's modulus of the 
pipe. 


Eq. (1) can now be written in the alternate form 

- pv < la > 


The following table gives v and -Ap/Av for some typical liquids in a 
rigid pipe ^ 



Velocity of Sound 
Y p , fps 

Pressure Increment. Ap, psi 
Velocity Reduction, -Av, fps 

Water (60°F) 

4700 

63.3 

Liquid Oxygen (-325°F) 

3370 

53.6 

Liquid Hydrogen (-428°F) 

4010 

3.92 


Within a pressure wave, the velocity of the fluid particle is given 
by Eq. (la) 


-Av = Ap/pv 


P 


It is in the direction of wave propagation if the pressure wave is a com- 
pression wave and is opposite to the direction of wave propagation if it 
is a rarefaction wave. 


The pressure wave will be reflected unchanged when it comes to the 
closed end of a pipeline. At the open end (where pressure is fixed) a 
pressure wave is reflected reversed in sign, i.e., a compression wave is 
reflected as a rarefaction wave and a rarefaction wave as a compression wave. 
At an obstruction (change in cross section of pipe, pipe bends, etc.), the 
pressure wave is partially transmitted and partially reflected. In particu- 
lar, it can be shown from Eq. (1) and the equation of continuity that at a 
cross section where the area changes from A^ to A an incident wave of 
intensity Ap is reflected as a wave of intensity TRef. 1) 



(4) 


A 1 - A 2 
A 1 + A 2 


Ap 


and transmitted as a wave of intensity 

2A. 


A 1 + A 2 


Ap 


With the above understanding we shall be able to discuss the 
of events following the opening of a valve. 
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(5) 

sequence 
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62.3 Instantaneous Valve Opening 


If a valve can be opened instantaneously, the fluid layer adjacent tfo 
the valve drops to the outside pressure with a resulting increase in velpcity 
Av related to the pressure drop -Ap by Eq. (1). This pressure drop, -Ap, 
propagates upstream with sonic velocity v • With the passage of this rare- 
faction wave, the fluid particles move downstream with velocity Av , At the 
upstream end of the pipe (a tank or pump outlet), the pressure is maintained 
constant, the incoming rarefaction wave -Ap is reflected as a compression 
wave of magnitude Ap moving downstream. The fluid particles now return to 
their previous pressure while acquiring a velocity 2Av. When this compres- 
sion wave reaches the valve end, it is again reflected as a rarefaction 
wave moving upstream. This process goes on and on until finally ^steady flow 
condition prevails. The following figure shows the idealized pressure and 
velocity variation at the mid-point of a pipeline. 


I 



Fig. 62,1 Pressure and Velocity Variation in a Pipeline 


The actual pressure- time and velocity diagram, of course, will have rounded 
corners and less stiff rises and falls because of friction and other effects. 
It is clear that in this case the pressure in the pipeline alternates 
between tank pressure and outside pressure, and there will be no dangerous 
high pressure build-up. This conclusion also holds if the opening of the 
valve is so rapid that the valve is fully open before the reflected compres- 
sion wave comes back to the valve, ice., 


Time of valve 


opening <= 


Twice the length of pipeline 
Velocity of pressure wave 
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62.4 Slower Valve Opening 


The preceding section indicates that a characteristic time of the pipe- 
line can be defined as 

T = Time for pressure wave to travel from the valve to the pipe 
inlet and back 

When the time of valve opening is less than T, we have the preceding case 
of instantaneous valve opening with no dangerous pressure build-up. When 
the time of valve opening is large compared with T, the elastic property of 
the fluid is negligible, and again there is no dangerous pressure build-up. 
In this section we shall discuss the case where the valve opening time is 
intermediate between these extremes and shall show that pressure surges can 
occur under certain (cases . 

In order to illustrate, we shall consider a valve partially opened 
instantaneously, with the opening held constant thereafter. The fluid layer 
behind the valve immediately experiences a pressure drop -Ap. The magnitude 
of Ap depends on percentage area of valve opening A. If the valve opening 
is not large, the usual discharge equation holds: 

Av = c d A"Y2p/p (6) 

where c^ is the discharge coefficient (close to unity) and 

P = P Q + Ap 

p being the pressure in the pipeline before valve opening. This equation 
tSgether with Eq. (la) 


-Ap = AvpVp 


(la) 


enables us to determine the pressure drop~Ap. This pressure drop “Ap (and 
velocity increment Av) propagates as a wave of rarefaction upstream and is 
reflected at the pipe inlet as a compression wave. After the passage ’of 
this reflected compression wave, the fluid returns to the pressure p and 
acquires a velocity 2v. When the reflected compression wave comes back to 
the partially opened valve, because of the obstruction of the valve, only 
part of the wave is transmitted, the rest is reflected as a compression 
wave of intensity 


Ap 



(7) 


according to Eq. (4). The pressure in the pipeline now jumps to the pressure 

+ ^(ts) (8) 

By keeping track of repeated reflections, we will be able to find out the 
pressure in the pipeline at each instant. The idealized pressure-time 
history at the valve end and at the mid-point along the pipeline is shown 
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in the following figures. The actual pressure- time history would, of course, 
be smoother because of friction effects and because of the impossibility of 
instantaneous valve opening,, 



I 

Valve opens time, t 


Fig.62 0 2a - Pressure-Time History at the Valve End 
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It is clear from these figures that each point along the pipeline experiences 
a decaying pressure oscillation of period 4i/v . The maximum jump in pres- 
sure 

i r - a i 

(9) 




1- A 
1+A 


is the same at points along the pipeline. But at points closer to the pipe 
inlet end, the duration of pressure quiescence is longer. Notice that since 
Ap^p Q and A^ 0, the pressure jump can never be greater than the static 
pressure of the fluid in the pipeline. 


As an example, we shall consider a pipeline carrying liquid oxygen at 
190 psig. A valve at 50 ft. from the pipe inlet is suddenly opened. The 
percentage area of opening is A = 1/4* For simplicity, we take the discharge 

coefficient c^ = 1 0 Eqs„ (la) and (6) become 


-Ap = 53.6Av 

AV , X 1 44 

from which we obtain the pressure drop following the valve opening 


(lb) 

(6a) 


-Ap = 189.1 psi 


The maximum pressure in the pipeline is 


P Q - Ap 


1-A 

1+A, 


= 190 + 189.1 


filial 


[1+1/ 4j 

The period of the pressure oscillation is 


303 psi 


44 _ 4 x 50 

v 33% 

P 


0.06 sec. 


In the above we considered the valve as being opened instantaneously 
and held constant thereafter. In the actual situation the valve is being 
opened continuously, and we have to know the variation of valve opening 
with time. The above example shows that most of the pressure drop occurs 
in the beginning of the valve opening; therefore, a good approximation td 
the maximum pressure jump can be obtained as 


P 


o 


1-A 


( 10 ) 


where X is the percentage area of valve opening at a time of — after the 
start of valve opening. v p 


A surge pressure of 325 psi was recorded by Marshall Space Flight Center 
following the opening of valve for liquid oxygen at 190 psi. This shows good 
agreement with the present result, even though the exact manner of valve 
opening is not known. Also recorded was the 50 milliseconds time interval 
between the peak and the valley of the pressure surges. This corresponds to 
thfe half period 2£/v of the pipeline. Assuming a sonic velocity of 
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Vp - 3370 fps, this would correspond roughly to a pipe length of 85 ft. 

It is to be noted that in the usual treatment of water hammer effect, 
the discharge equation (6) is used instead of Eq. (7) in the determination 
of the pressure surges. Water hammer is a transient effect. It seems to be 
unreasonable to use Eq. (6) which is based on the steady flow Bernoulli* s 
Equation. Furthermore, the use of Eq. (6) cannot explain the sequence of 
events following the instantaneous valve opening. It Is also in conflict 
with experimental data obtained in Huntsville and elsewhere. (The use of 
Eq. (6) can only predict pressure surges smaller than a quarter of the static 
pressure (Ref. 2),) The use of Eq. (7) in the present treatment seems to 
be more reasonable, based on physical reasoning as well as providing better 
correlation with experimental evidence. 
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62. 5 Closing of Valve 

The water hammer effect in the sudden closing of a valve is a subject 
treated extensively in the literature. We shall only describe it briefly 
here . 


When a valve is closed instantaneously, the layer of fluid in contact 
with the upstream side of the valve is suddenly stopped; its kinetic energy 
is converted into strain energy with a sudden rise in pressure. This gen- 
erates a compression wave moving upstream with sonic velocity, and with the 
passage of the compression wave, the fluid originally flowing with velocity v 
comes to a rest. Similarly, a rarefaction wave moves downstream of the valve. 
The amplitude of the compression wave can be obtained from Eq. (la) as 


Ap = -Avpv 


P 


and it can be seen that dangerous high pressures may occur. The propagation 
and reflection of these pressure waves follow much the same manner we dis- 
cussed previously. If the valve is not closed, instantaneously, the pressure 
in the pipeline will be the superposition of these traveling waves and will 
in general be lower than that given by Eq. (la). If, however, the valve 
closing time is shorter than the characteristic time T of the pipeline, then 
there is no chance for the reflected rarefaction wave to come back to the 
valve and prevent the pressure at the valve from reaching its maximum value 
given by Eq. (la). 
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62 . 6 Fluid Reaction at Pipe Bends 


We have shown that due to sudden initiation of flow, pressure 
surges approaching twice the operating pressure may occur. The increase in 
"hoop stress" in the pipe due to this pressure surge is usually considered 
in piping design. Because of the unavoidable existence of pipe bends, how - 
ever, large bending moment and axial force may occur, depending on how the 
piping is supported. There are two kinds of fluid reaction at a pipe bend, 
that due to the fluid momentum and that due to fluid pressure. The fluid 
reaction due to the change of fluid momentum can be obtained as follows: 



FIGURE 62.3 Reaction due to Change of Fluid 
Momentum at a Pipe Bend 


From Newton's law we have 

F = 2 Q v cos 9 

where Q = mass flow rate, slugs/sec. 

v = fluid velocity, ft. /sec. 

F = reactive force, lb. 

As an example consider a 90° bend in a 20 in. inside diameter pipe. The 
fluid is LOX at a pressure of 190 psi, a density of 2.3 slug/ft^ and at a 
velocity of 25 ft/sec. 


then F = 2 (2.3) (*)( -^-fy ) 2 (25) cos 45° 

= 180 lb. 

which is rather small in comparison with the end force of 

(190) <*)(—) 2 = 59,600 lb. 

due to static fluid pressure. This example shows that the fluid momentum 
effect is generally not of importance unless there is possible resonance due 
to pulsating flow. The static pressure on the other hand is of great im- 
portance for fluids under high pressure. Depending on the method of 
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supporting the piping, the static pressure force may be transmitted as an 
axial load, a bending moment on the pipe or may be transmitted entirely to 
the support. The axial load due to static pressure has been well recognized 
by previous investigators. The possibility of large bending moment, however, 
has not attracted much attention. The following figure illustrates the 
situation where the static pressure at a pipe. bend gives rise to bending 
moment but negligible axial load. 


expansion joint reaction due to static pressure, F 



FIGURE 62.4 Pipe Arrangement Giving Rise to Large Bending Moment 


Because of the flexibility of the expansion joint some of the hydraulic end 
load is transmitted to the support A. The pipe at A is subjected to a bend- 
ing moment of F Z. Take our previous example, and assume the unbalanced hydraulic end 
load is half the maximum possible with complete flexibility at the expansion 
j ointy 

F = 1/2 (59,600) lb. = 29,800 lb. 

a length Z = 2 ft., the bending moment at A is 
29,800 x 2(12) = 720,000 lb- lb 
which could be rather serious. 

So far we have been considering the hydraulic end load under steady 
flow conditions. With the sudden intiation of flow, we have shown that surge 
pressures higher than the operating pressure may occur. First. of all 
this higher pressure aggravates the adverse effects of hydraulic end load we have 
discussed above. Secondly it is a source of shock loading which may set up 
undesirable vibrations. The following figure illustrates a pipe arrangement 
which is subjected to negligible bending moment during steady operating 
conditions but to a shock loading during start-ups. 
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FIGURE 62.5 Pipe Arrangement Giving Rise to Impulse 
Following Valve Opening 


Because of the pressure oscillation in the pipeline after sudden opening of 
the valve, the pressure at A and B may not be the same at any instant. This sets 
up an unbalanced force in the pipe. To illustrate, consider the valve at C 
being opened instantaneously. The pressure drops from the static pressure of 
approximately 190 psi to, say, vacuum. This pressure drop propagates up- 
stream. After the wave front passes B, the hydraulic end force acting toward 
the left is reduced to zero, while at A the hydraulic end force acting toward 
the right is still 190 psi„ There is an unbalanced force of (190) (pipe 
cross-sectional area) = 59,600 lb. acting to the right, which is taken by 
the pipe support. After a time interval of 


_ pipe length between A and B 

velocity of wave propagation 

the wave front reaches A and the unbalanced force vanishes. Since the wave 
velocity is usually large, this time interval is small and this unbalanced 
force appears as a sharp impulse 

P 

59600 lb 

_L 

FIGURE 62.6 Impulse Due to Fluid Pressure 

which may excite undesirable vibrations of the pipe system. This pressure 
wave is reflected back and forth and, therefore, will appear repeatedly until 
it finally dissipates. The situation becomes most , serious when there is 
possible resonance . 
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62. 7 Conclusions 


In this section the pressure surges in a pipeline following the sudden 
opening of a valve have been investigated. It is concluded that 

(1) There will be pressure oscillations in the pipeline with period 
equal to 

k!L _ 4 (pipe length from valve to pipe inlet) 

v Velocity of wave propagation in the pipeline 

P 

If the valve opening is either very fast or very slow, there will 
be no high pressure build-up. For valve opening time intermediate 
between these extremes, there will be high pressure surges. The 
maximum amplitude will always be smaller than the static pressure 

p In the pipeline before opening and can be approximated by 
o 



where A is the percentage area .of valve opening at time 2^/v after 
the start of valve opening. At points closer to the pipe inlet, 
the duration of the high pressure is shorter. 

(2) These high pressure surges may affect the sealing capability of 
the fluid connect©* in the following ways: 

(a) Th^ aggravate the adverse effects of internal pressure such 
as stretching the pipe line and reducing the gasket com- 
pression, causing the flanges to roll due to the barreling 
effect, and introducing higher stresses in the pipeline. 

(b) Ttey ace a source of shock loading which may induce undesirable 
vibration of the pipeline. 

(3) Recommendations : 

(a) The usual proof pressure should be raised from 1.5 to 2 times 
the operating pressure, since, as shown in Eq. (10), the 
surge pressure may be greater than 1.5 times the operating 
pressure but never greater than twice the operating pressure. 

(b) It is desirable to locate the fluid connector as far as 
possible from the valves, since the duration of pressure 
surges will be shorter there. 

(c) Pipe layout should be treated with care to minimize the 
conditions that give rise to large bending moment due to 
hydraulic end load and the conditions that are prone to shock 
loading during start-tip and shut-down. 
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EFFECT OF SHOCK AND VIBRATION ON CONNECTORS 


by 

G. W. Sarney 


63.0 Summary 

This section considers a simplified analysis of the missile piping 
system to determine its response in terms of pipe reactions of interest to 
the fluid connector designer. The sinusoidal and shock excitation levels 
specified for the Saturn C-l, Block II vehicle are used on typical pipe con- 
figurations to predict their response in terms of end bending moments. A 
dynamical method of analysis is presented for simple pipe configurations to 
show the nature of the problem. For complex pipe systems a simplified 
approach is used to give the fluid connector designer a quick estimate of 
the pipe reactions. The loads determined in typical examples are quite high 
and must be included in the connector design if leakage is to be avoided. 
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63 . 1 Introduction 


The fluid connector joins sections of pipe and seals against leak- 
age. Environmental conditions in the missile induce vibrations in the piping 
system which are transmitted through the fluid connector (Fig. 63.1). Design 
of the fluid connector requires an estimation and evaluation of the effect of 
these pipe reactions. 


The random noise excitation may be predicted based on past tests 
and is presented in terms of vibration acceleration levels in the Shock and 
Vibration Specifications (Ref. 1). This loading acts on the pipe configura- 
tion which may be described by its end conditions, joints, physical dimensions, 
and damping characteristics. The response of the system can then be predicted 
with varying degrees of sophistication. This analysis uses a greatly simpli- 
fied approach to meet the needs of the fluid connector designer to obtain 
estimates of the missile piping end reactions. 
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63. 2 Environment 


The piping system of a missile is subjected to vibrations due to such 
environmental effects as maneuvering loads, aerodynamic noise, the rocket 
engine sound field, and direct transmission from the supporting structure. 

These effects are grouped and presented in the "Shock and Vibration Specifi- 
cations for Components on Saturn, C-l, Block II Vehicles" (Ref. 1) . This 
specification presents acceleration loading for sinusoidal resonant vibrations 
in terms of peak "g" level for a given ran§e ( of frequency. The shock test 
requirement specifies an impulse shape and g level of shock acting along any 
axis of the component. 
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SATURN BLOCK II, LOX DUCTING 




63.3 Definition of Component Loading 


The vibration excitation in the missile environment is random 
in direction, magnitude and frequency. This random excitation acts on a 
piping system which usually has many intermediate joints connecting different 
shaped sections (see Fig. 63.1). A complete dynamical analysis of such a 
system is not required for the purposes of the fluid connector designer. The 
procedure here, then, is to make reasonable simplifications of the loading 
and structure to obtain quick estimates of the pipe end reactions. 

The following simplifications of loading are made. The 
vibration acceleration level specified for the component under consideration 
will be the value for the frequency range which includes the natural 
frequency of the beam. This acceleration acts on the mass of the beam to 
produce a uniform load transverse £o the beam and a concentrated load along 
the beam axis. The most severe vibrations will occur in a direction which 
is transverse to the longest beam length. 

Complex pipe configurations will be resolved into straight beams 
whose length equals the pipe length normal (transverse) to the direction in 
which the acceleration acts. Sections of pipe along the acceleration axis will 
serve only to produce a concentrated load. 
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63, 4 Response 


The response of several typical pipe configurations will be ob- 
tained for the loading conditions discussed. 

A simple configuration of a straight pipe with fixed ends (the 
actual end support may vary from simple- support to fixed; however, the effect 
of this variation on connector bending "moment will depend on connector 
location) is first considered and the response found in terms of the normal 
modes of vibration of a simple beam based on the analysis of Sections 63.5 
and 63,8. This analysis is useful to show the type of response typical of a 
missile piping system. 

An example of a complex piping system is included to illustrate 
the simplifications of loading and pipe configuration used to estimate the 
response (Section 63.6). A static analysis is used on this type of system 
to compute the inertial reactions of the pipes for a given "g" level of ex- 
citation. The pipe reactions thus obtained are modified by an amplification 
factor for possible resonances. 

The piping structure will contribute an appreciable amount of 
damping. The gimbal is a sliding type joint which serves as an additional 
damper so that the expected amplification factor at resonance, Q, has the 
range 


5<Q<20 


The conservative estimate Q = 20 will be used in the examples that follow. 



63. 5 Vibration Analysis of Lox Interconnector Line 


This first example is a straight pipe to which a complete analysis 
can be simply applied. 

The LOX interconnector line connects the main oxygen tank with 
smaller oxygen tanks on the periphery of the Saturn Block II Vehicle. 

Configuration (See Fig. 63.2) 


74 




FIGURE 63.2 LOX INTERCONNECTOR LINE 









63.5.1 Physical Parameters 


Internal pressure 
Density LOX 
Density Steel 
Beam Thickness 
Beam Diameter 
Beam Length 
Cross Section, Inertia 
Modulus of Elasticity 
Beam Mass per Unit 
Length 

Beam Weight per Unit 
Length 


200 psig (maxi 
p L - 62.4 lb/ft J 
p s = .285 lb/in 3 
t - .0625 in 
D q = 12.75 in „ 
i = 58.3 in: i = 3390 in 
I = 43.2 in 4 
E = 30,000,000 lb /in 2 

p = .0132 lb sec 2 /in 2 


U = 5.09 lb/in 


2 . 


11 . 5xl0 6 in 4 
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63.5.2 Excitation 


The LOX interconnector line is located in Zone 3-2-B as designated 
by the "Vibration and Shock Specifications for Components on Saturn, C-l 
Block II Vehicles" (Ref. 1) . In this zone the following vibration and shock 
requirements must be met along any axis of the IX)X interconnector component. 

Resonant Test Code 1-A 

20-28 cps at 0.1 in D.A. displacement 
28-72 cps at 4.0 g peak 
72-140 cps at 0.015 in D.A. displacement 
140-2000 cps at 15,0 g peak 

Shock Test Code b 

50 g peak 
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63 . 5.3 


Analysis of Loading 

All loading cases are considered with LOX in the pipes since this 
lowers the natural frequency and increases the magnitude of the load signifi 
cantly . 

The beam may be analyzed in either of the two models shown below, 
depending on stiffness of the gimbal joint. 





63 5 4 Response for Continuous Model 

The following mode parameters are obtained from Ref. 2 for the 
continuous beam with clamped-clamped supports. 


Mode 

3 n_ 

V 

(v ) 2 

(v ) 4 

f 

n 

0 n (x=£/2) 

0 n (x=O) 

n 

1 

0.98 

4.73 

22.4 

500.6 

330cps 

1.588 

2.00 

3 

0.99 

10.99 

120.9 

14617 

1790cps 

-1.406 

2.00 

5 

0.99 

17.28 

298.6 

89135 

4400cps 

1.414 

2.00 


63.5.4.1 Steady State Vibration 

The steady state acceleration loading in the resonant test code 
requirement is most critical in the frequency range of 140 to 2000 cps where 
the peak acceleration is G = 15 g's. 

This acceleration acts on the mass of the beam to produce the load 
intensity P(x) . (See Appendix, Sec. 63.8) 

P(x) = pG 


Then 


P = jj 
n & 


(pG)0 n (x)dx 


The integral is evaluated using Ref. 3, 


clamped-clamped beam. Integral #1 . 


P 


n 


(pG) 2a n 

~Tb 

n 


i-(-D n 


n = 1,2,3 


Since a j-,-1 and for n = 1,3,5, 

P = 
n^Bn 


We calculate the bending moment at the end of the beam (x-0) according to 
the relation obtained in the Appendix, Sec. 63.8 


M 


oo 



i 2p A (i!) 




n = 1,3,5 

(even modes are 
not excited) 


(72) 
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The amplitudes of the fifth and higher modes are negligible and their natural 
frequencies are too high for resonances to occur in the range of interest. 

The amplification factor is taken to be Q=20 . 

The amplitude of the bending moment is: 



at f = 330 cps we get the maximum moment ignoring phase relations as: 


M = 392,000 + 1600 = 393,600 in lb 


63,5.4.2 Shock Response for the Continuous Model 


The transient excitation requirements in the shock test code are 
most critical for the rectangular type pulse. The ratio of response to 
excitation which is applied to each mode of beam vibration is: (See 

Appendix, Sec. 6 3.8.3 


Response 


Excitation 


-fit 

e cos 


u t 


n 


0<t<«? 


for rectangular pulse T =.005 seconds 
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A calculation is now made to determine the time needed to reach the 
maximum response amplitude for each mode in this beam. The maximum amplitude 
of the shock response occurs when 

COSO) t “ -1 

n 

Then the rise time is: t = Jt/u) n 


Mode Circular Frequency, Rise-Time, t R 


1 (330) (6.28) rad/sec 

3 (1790) (6.28) 

5 (4400) (6.28) 


.0015 second 

.00028 

.00011 


Since there is little damping of the response during the excitation period 
( V = 5 milliseconds) it is possible for the first fives modes to reach 
their maximum during this time and this summation must be considered. 


The load function for a G » 50 g's peak impulse is 
g 





Now the load per unit length is: 


P 

n 


-if 

jo 


Si 


(pG )0 (x) dx 
n 


Using the integral tables of Ref. 3 for the clamped-clamped beam: 


P 

n 


(HG) (2a n ) 

_T \ 



n = 1, 2,3,4- 


Since a ; pG = (5.09) (50) 

1018 


n = 1,3,5 
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The bending moment at the ends of the LOX pipe for a 50 g shock impulse as 
determined in Section 63.8 is: 


00 t P 0 (x) 
M- I 


n=l 




-r| t 
e cosu) 


'-‘J 


The maximum amplitude of the end bending moment is therefore 
2 " 

<ne Z P 0 (x) j- -J 

M = Y — (2) n = 1,3 (even modes are not 

n=l (£B n ) 2 L J excited) 


Then the value of bending moment neglecting phase is: (higher modes are 

negligible) 

M = 131,000 + 10,400 = 141,400 in lb 
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63 5 5 Response of Discrete Model 


In this section the LOX interconnector is analyzed by considering 
the gimbal joint to have zero stiffness in bending. 



The beam loading is analyzed in two sections; however, the natural frequency 
of the system is a function of the deflections of the combined system. Since 
Zh, is 12 inches, the split in the continuous beam ( J ■ 58 inches) occurs 
at its point of inflection or when the moment of the continuous beam is zero. 
This discrete model will have the same static deflection shape as the con- 
tinuous model and hence nearly the same natural frequencies. 

An approximate type of analysis will be used. The cantilver beam 
Z^ will be analyzed as loaded with a uniform load and a concentrated load 
at the end due to one-half the loading on the center span The mode 

parameters will be taken for the cantilever section alone as obtained from 
Ref. 2 for the uniform clamped- free beam. These are as follows: 


Mode 

f * 
n 

*n<\> 

0/(0) 

fi* JB 
1 n 

<M n ) 2 

a n 

i 

330 

2.00 

2,00 

1.875 

3.516 

.734 

2 

910 

-2.00 

2.00 

4.694 

22.03 

1.018 

3 

1790 

2.00 

2.00 

7.855 

61.70 

.999 

4 

2950 

-2.00 

2.00 

10.995 

120.90 

1.000 

5 

4400 

2.00 

2.00 

14.137 

199.86 

.999 


* f for continuous beam model 
n 


63.5.5.1 Resonant Vibration 

Steady-state vibration response for a sinusoidal excitation of 
G = 15 g's is now calculated: 

Load: W = 1/2(5.09) (58-24) (G) = 86.6 G 

then W = 1300 lb. 

Then the load intensity is: 
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P(x) = W(x=/5 1 ) + |iG 


and 


W 


P n -ir«„ <i l> + 


IJ.G0 (x) dx 
• o n 


Using the tables of integrals in Ref. 3 for a cantilever beam 


p „ - \ W +2 


*l B n 


n = 1, 2, 3,4,5,- 


The steady state response for moments at the fixed end of the cantilever 
beam based on the relation from Section 63.8 is 



Neglecting phase, as a conservative assumption, the amplitude of this bending 
moment for the first four modes then is 






The amplification factor at resonance is taken to be Q— 20 as for the continuous 
beam: 
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The amplitude of the end bending moment at the worst condition when, f = 330 cps, 
is: 


M - 454,000 + 2,000 + 1,000 = 457,000 in lb. 


63.5.5,2 Transient Solution for Cantilever Beam 

The analysis of the response function for a rectangular shock pulse 
is similar to the approach used for the continuous model. Therefore, the 
maximum ratio of response to excitation is two for each mode. The excitation 
function for a G = 50 g's rectangular impulse on the cantilever beam is 
similar in form to the steady state vibration excitation. 


P = 
n 


MOL 


tt W +2 n 


1 n 


The maximum amplitude of the clamped end bending moment therefore is: 


M 


Z * P 0 (x=0) 

_0- I (2) n * ■ — 
n=1 


Considering the first three modes, and neglecting phase: 


M = 151,000 + 16,000 + 7,300 = 174,000 in lb. 
x=0 
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63,6 


Vibration Analysis of LOX Suction Line 


The LOX suction line is an example of a complex pipe configuration 
which must be simplified. 

The LOX suction line connects the LOX tanks to the nozzle assembly. 
A typical configuration from the Saturn Block II Vehicle (Fig. 63.1) is shown 
in the attached drawing (Fig. 63.3). 
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LOX SUCTION LINE - TOP VIEW 

FIGURE 63.3(a) 
63-20 
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LOX SUCTION LINE - SIDE VIEW 
FIGURE 63.3 (b) 

63 21 


88 



63.6.1 Physical Parameters 

The following parameters 
continuous in the system: 

Density LOX 

Density steel p g 

Line thickness t 

Line diameter D q 

Cross section inertia I 

Modulus of Elasticity E 

Line mass per unit 

length p 

Line weight per unit 

length p 


are typical and are applied as though 

= 62.4 lb/ft 3 

= .285 lb /in 3 

= .0625 in 

= 8.00 

= 11.8 in 4 

= 30xl0 6 psi 

2 2 

= .00573 lb sec /in 
= 2.21 lb/in 
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63.6.2 Loading Requirement 


The LOX suction line is located in Zone 1-1-G as designated by the 
"Vibration and Shock Specifications for Components on Saturn, C-I, Block II 
Vehicles". In this zone the following steady-state vibration and transient 
shock conditions must be met when the inputs are applied to the end of the 
line which fastens to the pump on the lower end. 


Resonant Test Code 1-A 

20-100 cps at 4.0 g peak 
100-142 cps at 0.0078 in D.A. displacement 
142-540 cps at 0.00055 in D.A. displacement 
540-1000 cps at 0.00055 in D.A. displacement 
1000-2000 cps at 28.0 g peak 

Shock Test Code b 

100 g peak 
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63.6.3 Analysis of Loaditi; 


The procedure here is to consider the system to have uniform stiff 
ness and to consider the acceleration to act along each of the three axes 
shown on the drawing. The pipe system is split into four sections labeled 
A, B, C and D which are either normal to or lie along the axes x, y and z. 
Each section has a static weight per unit length given by p = 2.21 lb/in. 

The complex pipe configuration will be simplified by constructing a beam of 
length equal to the pipe length which is normal to the direction of the 
acceleration. The loading diagram for acceleration acting along the y axis 


(Section C) W - (2 . 21) (42. 2 )G 



The pipes which are normal in direction to the axis through which the 
acceleration acts contribute a uniform load along their length, and pipes 
which are parallel to the acceleration axis contribute a concentrated load. 
The bending moments at the fixed ends can be calculated for a steady load 
distribution at the "g" level of excitation along each of the three axes. 
These values are modified by an amplification factor for possible resonances 
(Q=20) . 

This procedure is used to find the bending moments at each end of 
the LOX suction line for the three cases: in which the acceleration is 
resolved along the x, y and z axes respectively. The accuracy depends on 
the ability to obtain the same deflection shape in the equivalent model with 
the new length and loading conditions. 
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63.6.4 Res P° nse 


The first case is analyzed using the loading diagram above for 
acceleration resolved along the y axis. The moments at the ends are evalu 
ated by superimposing the following two loading diagrams and applying 
formulas from Roark, Ref. 5. 


Then 



M 


2 


' 2.21(77)! 
12 


+ (62.7) (14.3) (42.2) (2.21) 

(77; 


G = 2020 G 


For this component G = 4 g's, assuming the natural frequency below 100 cps 


M x = 5,200 

in 

lb 

M 2 = 8,000 

in 

lb 


Applying the expected amplification factor for resonance 
M 1 = 104,000 in lb 
M 2 = 161,000 in lb 
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The second case is analyzed using the following loading diagram 
for weight resolved along the x axis. 

(Section B) W = (33)(2.21)G 





For the LOX suction line G = 4 g r s 
M 1 = 6,800 in lb 
M 2 = 6,000 in lb 
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The end pipe reactions at resonance are then: (Q * 20) 


M x - 136,000 in lb 
M 2 - 120,000 in lb 

The third case is analyzed using the following loading diagram for weight 
resolved along the z axis. 



Nozzle 

End 


The static bending moments at each end are: 

M x = ( 2 - 21 ) ( 89 -5 I G = 1550 G 
M 2 = 1550 G 

For the LOX suction line G = 4 g's 


= 6,200 

in 

lb 

M 2 = 6,200 

in 

lb 


Applying the amplification factor for resonance } Q=20: 
= 124,000 in lb 
M 2 = 124,000 in lb 
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63,7 Conclusions 


Only for simple pipe configurations can an exact dynamical 
analysis be easily applied. For pipe systems with intermediate joints or 
curved sections, the normal modes of vibration cannot be predicted easily, 
and simplifications must be made. A pipe system of variable stiffness, com- 
plex shape and uniform loading is replaced by an equivalent straight beam 
with a uniform stiffness, an altered length and a variable loading to ap- 
proximate the deflection curve of the original configuration. The typical 
examples of this section show that gimbal-type joints of indeterminate 
stiffness occur near the point of in. action of the continuous equivalent 
system so that a good approximation for the deflection curve of the original 
configuration is obtained. A static analysis of the equivalent system can 
be used to obtain a reasonable estimate of the pipe end reactions for 
purposes of the fluid-connector designer. 

The end reactions determined on some typical pipe configurations 
in the environment of the Saturn C- 1 Vehicle are quite high and must be 
considered in the connector design if leakage is to be avoided. 
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Appendix: Vibration Analysis of Beam 


3 , 8,1 Normal-Mode Analysis (Summary of References 3 and 7) 

The governing differential equation of motion of a beam with a 
time-varying load is obtained by equating the external load intensity to 
the sum of the inertial and elastic loads. 


Y 



The equation of motion is written with respect to the co-ordinate system, X,Y 
Motion of the beam with respect to this co-ordinate system denoted by "y" may 
be caused by an application of an external load or by movement of the beam 
supports which is denoted by Y^ and X^- The movement of the supports with 
respect to the fixed co-ordinate system produces an equivalent "G" loading on 
the beam. This "G" loading is a continuous load which has a value at x=0 of 

A x 

P(0) =p — o' CD 

St 2 


and at x 

p(^) 



( 2 ) 


This equivalent load intensity varies linearly across the beam to 

satisfy 

the boundary values. 9 

^ 2 


a y x 

a z Y 


For the case of uniform acceleration _ - 

2 - Gg 

(3) 

at 

at 2 


the equivalent load intensity becomes 

p (x, t) = pGgf(t) 


(4) 

or 

P (x, t) = |iGf (t) 


(5) 
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where 


= mass per unit length of beam 
p = weight per unit length 

G = number of g's - uniform acceleration 

f(*-) = time function of acceleration loading 

The external load produces a beam deflection y with respect to the ends of 
the beam. This deflection is resisted by the following forces. 

1, Flexural stiffness of beam - resistance to bending along beam. 


S y 

M - - El 

> 

bending moment 

(6) 

Sx 2 




Q " S * 


shear force 

(7) 

FC*,t> = 

El 

EI A > 

load intensity 

(8) 

Sx 4 




where: El is uniform stiffness of beam 

2. Beam Inertia - resistance to movement of each section 

P(x,t) = p (9) 

St 2 

where: p = uniform density of beam 

3. Damping in beam - considered viscous and distributed along the beam 

P(*,t) - c §f (10) 


The governing equation for 
shear forces and rotary inertia, is: 


flexural vibration of the beam, neglecting 


El 


Sx" 


S 

St 


+ p “ + C St 


_ 


pGf(t) 


where : 


y = deflection of beam with respect to its ends, (in) 
El = bending stiffness of beam - uniform (lb/in ) (in^) 

p = uniform mass density per unit length (lb/sec/in ) 

C = damping coefficient per unit length 
p = beam weight per unit length 

G = number of "g's" peak acceleration 

f(t)= time function of acceleration loading 


( 11 ) 


( 12 ) 


The beam is a continuous system with an infinite number of degrees 
of freedom which may be analyzed by superimposing the response in the normal 
modes . 
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Let beam deflection in each mode y = 0 (x) 7 (t) (13) 

n n 'n 

Oo 

Then beam deflection y = ^ 0^ (x) 7 (t) (14) 

n=l 

where the mode shape in the nth mode is 0 ^ (x) and the amplitude 7 ^. 


The normal modes are solutions of: 



Trying as an assumed solution 

y = 0 (x) sino) t 
J n n n 

then Eq, (15) yields: 


Define 


then 


Define 


then 


Sx 


«n (x) 


p 03 

1 ~E~ K <*> ■ 0 


n 


P 


n_ 

El 



0 (x) - 0 = 0 

n n 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 
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The general solution for the mode shape 0^(x) from equation (19) is: 

0 (x) = AsinB x+Bcosp x+Csinhp x+Dcoshp x (22) 

n n n n n 

The constants A, B, C, and D are solved from the boundary conditions of mode 

shape at the ends of the beam and may be grouped as follows: 


0 n (x) - coshp n x-cosP n x-a n (sinhp n x-sinP n x) (23) 

The mode shape for the case of a simply supported beam will be analyzed here 
to determine p n , and 0 ( x ) . These values are tabulated for all other 
mounting conditions for tne first five modes in Ref, 2. 


Boundary conditions - simply supported beam 

y = 0 ; x = 0, x ® i (24) 

y n = 0; x = 0, x * i (25) 

Applying the above conditions to Eq. (23) yields: 

0 (x) * sinp x (26) 

n 7 n 

Using y = 0; x = i (27) 

0 - sinp Jt (28) 

r n 

then 

B Z = nit n - 1,2,3, (29) 

r n 


From Equation (18) 



GJ n = (nW-^S/El/p 1 

(30) 

and 

0 n (x) = sin mtx/^ 

(31) 


The response of the beam system will now be obtained for a transient 
and a steady-state type time-varying load. 
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63,8.2 Response to Steady-State Vibration 


The steady-state excitation is a sinusoidal force which is specified 
by its frequency and peak acceleration in g's. 


Excitation P(x)sincot * P 0 (x) sincot 

n n 

n=l 

r£ 

where: P = — f P(x)0 (x)dx 


n i 


n 


Response 


y = 


oo 

/ a 0 (x) sinuit 
* — n n 


(32) 


(33) 


(34) 


n=l 


Substituting the excitation and response into the governing differential 
equation yields for each mode n. 


EI 3n "<* w 2 a n 0 n (x)+Jfa,Ca n a n (x) = Vn (x) (35) 


where 


then 


s7p» (x) ] ’ 7 ®" <x) 

Ox L J P 


2 4 

CO = EI p /p 
n n 


2 2 

a co - a co + j ccoa /p * P / p 
nn n n ^ n M 


n 


3 n . 2 2. 

p(co^-co )+JC0C 


and 


<30 P 0 (x) sin(cot-V) 
*r n n 

n=l 


2 , 2 2 * 2 2 

p (co -co ) +c CO 


where tan ^ 


2 2 

= coc/p (co^-co ) 


(19) 


(18) 


(36) 

(37) 


( 38 ) 
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( 18 ) 


2 4 

Making the substitution co^pF Eip^ 


= £ 

rt=l. 


P (t (x) sin(tot- 'I) 
n n 


(39) 


Eip 



which can be written as 

oo 


£ Vn22.il 

n=l cs n ) 4 EI 


sin(<jot-jO 



A 2 


1 - ^1 + 


co 



CO 


300 


CO 


n * ^ n 


The deflection is the infinite sum of the modal deflections due to an excitation 


Z 

P n - ff P(x)0 n (x)dx 

acting on a mode shape characterized by 0^(x) an ^ P n * Integrals of this 
type are tabulated in Ref. 3. 

The beam response is also a function of the beam parameters: 

Z = length 
EI = stiffness 
p = density 

C = damping coefficient 

The amplification factor of a particular mode of resonance, co = co^, is of 
interest . 


Resonant deflection y = Q y L . 

•'n static 

Resonant condition co = co 

n 

Then combining equations (40) , (41) , and (42) 



The damping factor may be defined as 

■n = C/2 p 


(41) 

(42) 


(43) 


(44) 
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Amplification Factor for Shock 


The shock response of the linear, mul tiple-degree-of-f reedom beam 
can be obtained by applying the transient response for a single-degree-of- 
freedom system to the normal coordinates of the beam system. The normal 
coordinates are the amplitudes in the modes. The amplitude of the steady- 
state deflection for each mode is: 

A = F P - — -ii- (45) 

n ) 4 El 

This amplitude of displacement of the beam for each mode "n" will be modified 
by the shock response. 

The shock requirement used on the LOX pipes is presented in 
'Vibration and Shock Specifications for Components on Saturn, C-l, Block II 
Vehicles" (Ref. 1) in terms of peak acceleration for each component or block 
within the vehicle. The shock test code states that the acceleration shall 
be applied along each of the three perpendicular axes and shall be one of 
the following wave shapes. 

1. Triangular pulse with duration of 10 milliseconds 

2. Half-sine pulse with duration of 8 milliseconds 

3. Rectangular pulse with duration of 5 milliseconds. 

This analysis considers lateral vibrations which are the most 
critical. The shock pulses have approximately equal areas and approximately 
equal responses. The response of a linear, single-degree-of-freedom system 
with damping for the rectangular shock pulse is now derived. (Ref. 6, pages 
8-49,50). 


Equation of motion 


mv+cv-fkv = g(t) 

(46) 

fv-lfv+ v = |(t) 

(47) 

!(t) = ^ 

(48) 


v = response 
t = excitation 
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The Laplace transform with respect to "t" is: 


-^[s 2 F v (s)-Sv(o)-v(o)] + | [sF v (s)-v(o)] + F v (s) = F^ (s) 

Cd V 

n 


The initial conditions in this case are: 
v(o) = 0 
v (o) = 0 


The subsidiary equation becomes 

F (s) 

F.,(s) ' — 2 * 


Vt S + l 

CO 

n 


Response 


v(t) = L 


-1 




s24 {^n) S + W n 


For rectangular pulse I (t) = £ , 0<^t^- f 


Pulse 



-*• t 


F^ (s) = L je (t). 


= t/S, 0< t < T 


Then 


/.X - 2 T- 1 

v(t) = | u n L 


i s (s 2+ fal ) s + 


(50) 

(51) 

(52) 

(53) 

(54) 

(55) 


(56) 


(57) 
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Roots of denominators, in addition to S = 0, are 


! cco 


S = - w 


n + 


n 2k 


2 2 
c w 

n 

4k 2 


- 1 


Let 


CO) 


2k 2roto 


Rewriting Eq. (57) : 


v(t) - (g / n )L l 


( C +i^ i “ 2 J] s -to n C5 - i7 1 -i 2 

Therefore from inverse Laplace tables 


. * r 2 
v(t) = lu n 


, -at -bt 
be -ae 


ab ab(a-b) 


where : 


w n (C -i/i-.;) 


- co n (^+i J 


i -/) 


then 



Therefore Eq. 

b (t) 


(61) becomes 
l+e” W n 





— sin(w J171 

^ “ 


( 58 ) 


(59) 


(60) 


(61) 
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Now 


(59) 


Let Ti = — damping factor 

4m 


then tj = £ u) 


Since the damping £ is small 


i 5 C 


and the frequency 


“d = a> n ly 1- C 


v(t) = I 1-e 1 (cosu^t) tl t 


Therefore the dimensionless amplification factor for shock in the "n” mode 
of vibration of the continuous beam is: 


7 n (t) = = 1 - e _rit (coso) n t) 


Then for 0 < t < ? the transient response is given by Eqs. (14, (45) and (67) as, 


0° P 0 (x) or 

n n * t .-It 
v = i- e cos co t 

n=1 (4J EI L 


lOo 
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The solution for bending moment in the vibrating beam based on the 
deflection relation follows: 

oo 

y * I A n (X) "n (t) (69) 

n=l 


where A n y (t) = time response of nth mode 

A 

Bending Moment M = -El — ~ 

dx 


M 




7 n (t) 


where 



II 

0^(x) Tabulated values in Ref. 2. 


Then 


Jt P 0 (x) 

h- nn 
n=1 


7 n (t) 


Solution for flexural stress in vibrating beam 

_ m 

_ey max 21 


{ 


M 

D 

I 


bending moment, in- lb 
diameter of pipe 

moment of inertia about neutral axis 
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(70) 


(71) 


(72) 


(73) 



Solution for Shear force in vibrating beam 


where 


then 


Shear force 



( 74 ) 



! t r 

0 n <x> 



9 n M 


3 


P 


3 

n 


i I r 

0 (x) Tabulated in Ref. 2 

n v 


Q - 


oo 

n= 1 


<*„> 


* n (t) 


(73) 


(76) 


10^ 


■k 


63-40 



63.9 


References 


1. "Vibration and Shock Specifications for Components on Saturn, 

C-l, Block II Vehicles," Internal Note by Dynamics Section, 

George C. Marshall Space Flight Center; IN-P VE-S-62-7; May 1962. 

2 . R. P. Felgar and D. Young, "Tables of Characteristic Functions 
Representing Normal Modes of Vibration of a Beam," The University 
of Texas Publication, Circular 44, July 1949. 

3. R. P. Felgar, "Formulas for Integrals Containing Characteristics 
Functions of a Vibrating Beam," The University of Texas Pub- 
lication; Circular No. 14, 1950. 

4. S. H. Crandall, Random Vibration , New York, John Wiley and Sons, 
Inc . , 1958. 

5. R. J. Roark, Formulas for Stress and Strain , New York, McGraw-Hill 
Book Company, Inc., 1954. 

6. C. M. Harris and C. E. Crede, Shock and Vibration Handbook , 
McGraw-Hill Book Company, Inc., 1961. 

7. S. Timoshenko, Vibration Problems in Engineering , D. Van Nostrand 
Company, Inc., 1937. 



63-41 


64. EFFECTS OF RADIATION, VACUUM, AND THERMAL ENVIRONMENTS 
ON POLYMERIC GASKET MATERIALS 


by 

J.R. McLoughlin 


64.0 Summary 

The properties of polymers which are important for their consideration 
as elements of seals are discussed. The effects of radiation, high vacuum, 
high temperature, and low temperature on the mechanical behavior of various 
types of polymers are described. Thermal expansion, the glass transition, 
stress relaxation, oxidation, thermal degradation, and dynamic behavior are 
considered. 


1.09 
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64. 1 Introduction 


The chief interest in polymers for use in seals is as the gasket 
material. In this application the most important properties to be considered 
are the mechanical properties of the polymer although resistance to attack 
and degradation by other materials in the environment is, of course, a 
necessary requirement. 

Ideally, a gasket material should respond instantaneously to applied 
pressure and exert back pressure against the flange or compression element 
indefinitely. Furthermore, the material should be able to conform readily 
to temperature changes in the seal area without releasing its back pressure 
or moving laterally with respect to the compression elements. Such ideal 
requirements are, of course, very difficult to achieve with real materials 
except over limited temperature regions. For example, a network polymer 
in that region of temperature where it exerts ideal rubber behavior will 
very nearly conform to these requirements. However, if the temperature is 
raised very far, the rubber will relax its back pressure due to degradation 
of the polymer network by oxidation or simple thermal breakdown. Even the 
most thermally resistant rubbgrs such as silicone rubbers will be attacked 
slowly at temperatures of 250 C, and the degradation proceeds more rapidly 
as the temperature is raised. This causes the stress in the rubber to relax 
so that eventually the gasket ceases to press back against the seal flanges 
unless the latter are designed to continuously produce pressure by spring 
action. This process which takes place in a rubber is known as chemical 
scission, and the stress relaxation which it produces is called chemical stress 
relaxation. Some rubbers undergo such behavior at lower temperatures than 
others. For example, polysulfide rubbers undergo chemical scission by a 
process of opening sulfur bonds at 60-100°C. Ordinary carbon-based natural 
and synthetic rubbers are somewhat better, depending on exact formulations and 
anti-oxidants used, but none are very durable at 150°C. More thermally 
resistant polymers are being investigated, but it is doubtful whether the 
temperature limit will be raised much above 300°C in the near future. Therefore, 
polymers can be considered useful as gasket materials only at temperatures of 
300 degrees Centigrade and lower. 

As the temperature is decreased, another characteristic behavior of 
polymers must be considered. This is the transition from a rubber-like 
material to a glass -like material. All polymers exhibit this behavior, 
although the characteristic hardening takes place at different temperatures 
for different materials, and the shape of the modulus -temperature curve may 
be somewhat different due to variations in crystallinity, polarity, chain 
length and other structural variations. This transition is known as the 
glass -transition. Most polymers exhibit this behavior at room temperature 
and above, but some hydrocarbon rubbers have transitions at temperatures well 
below 0°C, such as polyisobutylene at -74°C, natural rubber at -73°C, and 
GR-S at -61°C. Some silicone rubbers also have very low glass transitions and 
are more suitable for use as rubbers at such temperatures than carbon-based 
rubbers because they do not have as great a tendency to crystallize. 

Below the glass transition, all polymers behave like hard solids with 
Young's moduli of elasticity in the neighborhood of 10 10 dynes/cm. 2 (145,000 
pounds/sq. inch) and higher. Above the glass transition, a polymer deforms 
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by uncoiling of t long chains of atoms which are linked together in a network 
of periodic lf crosslinks u . This is what makes it possible to stretch a rubber 
several hundred percent. The modulus may be as low as 100 P 0 S,I, As the 
rubber is cooled through the transition, it is no longer possible for the 
chain segments to move freely, and finally they cannot move at all, so that 
any deformation due to stress can only come through small displacements of 
chain atoms with respect to neighboring atoms on the same or nearby chains. 
Because this "freezing- in" process occurs over a rather wide temperature in- 
terval, there is another type of stress relaxation which occurs in high poly- 
mers at and below the glass transition (and sometimes in crystal transforma- 
tions). This relaxation is associated with short movements of segments of 
chains, and its chief importance for this study is that it results in a slow 
decay of the reaction force which is exerted by a gasket polymer on its re- 
straining flanges. Considerable data on this stress relaxation of polymers 
in their glassy state (below the glass transition temperature) have appeared 
in the literature, but most of the direct measurements of relaxation do not 
extend far below the glass transition, due to the small deformations and long 
times which are involved in such studies. Instead, the majority of investi- 
gators have used dynamic measurements (vibrations) to study the deformation 
behavior of polymers at low temperatures. Such measurements have been ex- 
tended down to as low as 4° Kelvin in some cases and to 80 Kelvin in the 
case, of quite a few polymers. A good review of the work done in this field 
was published recently by A 0 E 0 Woodward and J,A. Sauer (Ref. 1), This article 
presents dynamic shear modulus data and dynamic loss data for such polymers as 
polymethacrylates, polyvinyl ethers, polyvinyl esters, polyethylene, poly- 
amides, polyurethanes, polyesters, polyvinyl chloride, polyvinylidene chloride, 
polytrif luorochloroethylene (Kel-F) , poly tetraf luoroe thylene (Teflon) and 
polyacrylonitrile , 


To understand the importance of these measurements to this project, it 
is only necessary to realize that the small- de flection behavior of a polymer 
in various modes of distortion such as creep, stress relaxation, and constant 
rate of strain can be estimated from its behavior in such dynamic measurements. 
A wealth of information has been published on this subject, much of which has 
been collected in J.D. Ferry’s book, Viscoelastic Properties of Polymers 
(Ref. 2). If we are interested in the probability of stress relaxation of a 
polymeric gasket at low temperature which may result in failure by excessive 
decrease in the force the gasket exerts against its restraining members, we 
must convert the dynamic data to s tress- relaxation data. This can be done 
approximately by looking at the dynamic loss at a given frequency and temp- 
erature, A peak in the dynamic loss- temperature curve means that a mode of 
relaxation occurs at that temperature. Conversely, a low value of dynamic 
loss indicates low relaxation at that temperature for that frequency. Polymer 
properties are very sensitive to time as well as temperature, so behavior at 
several frequencies is necessary to gain a true understanding. However, it is 
usually possible to make a frequency- temperature conversion so that we can 
predict low-frequency behavior from high-frequency behavior by sliding the 
curves along the temperature scale. That is, the same behavior observed at a 
certain frequency and temperature will be observed at a lower frequency at a 
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lower temperature. Considerable experience is available for making this 
time- temperature transformation based on many polymers. So, from the pub- 
lished information on dynamic losses of polymers at low temperatures, there 
is available a body of information enabling us to make an approximate pre- 
diction of long-time stress relaxation behavior at low temperatures. In 
general, this shows that most polymers have very little dynamic loss below 
around -100 C. The most notable exceptions are those polymers which appar- 
ently exhibit some crystal transformations at lower temperature, such as 
polyethylene, polytetraf luoroethy lene (Teflon), and polyamides (Nylon), Even 
these have disappeared at -200 C. Therefore, we can expect very little stress 
relaxation in polymers at temperatures of -200 C 0 Instead they will behave 
like hard solids with a high modulus of elasticity and typically will tolerate 
only low strains without failure, sometimes less than 1 %„ 

One of the most serious limitations in the use of polymers as gaskets 
lies in the high coefficient of thermal expansion of polymers compared to 
metals which are likely to be used in the restraining flanges and compression 
elements. Polymers have a higher coefficient of expansion in the rubbery state, 
above the glass transition, where values of linear expansion around 1 to 2 
x 10 per °C are normal. This is about 10 to 20 times as high as most metals. 
The thermal expansion coefficients of plastics are only about 1/2 to 1/3 as 
high below the glass transition as above it, but they are still much higher 
than for metals. The discrepancy between metals and rubbers is not so 
serious, s ince a rubbery gasket can readily distort to conform to its metal 
retaining element. However, as it is cooled below the glass transition, the 
capacity for the plastic to deform and match the metal surface is diminished 
rapidly and is practically non-existent at very low temperature. Assuming 
a seal held in place by metal compressing elements on top and bottom surfaces 
and cooled to low temperature after compression, the reactive force by the 
plastic against the compression elements will begin to decrease as soon as the 
temperature drops through the glass transition. It may actually decrease to 
zero and the gasket may shrink away from the compressing faces if some provision 
is not made for a "follow up" spring action by these compression faces. While 
the flanges can be designed to achieve such a spring action in the direction 
normal to the compression surface, the mismatch of expansion in a direction 
parallel to the compression surface is much more difficult to compensate, since 
the gasxet and compression face are in constant contact. This may result in 
sliding of the gasket across the compression face so that irregularities in 
the face no longer match the irregularities in the plastic gasket which were 
molded in at much higher temperature. 


These considerations suggest that if relative sliding motion proves to be 
a problem, the coefficients of expansion of compression face and gasket might 
be matched by incorporating mineral fillers in the polymer to lower its 
coefficient of expansion or by making the retaining mechanism out of reinforced 
plastic in order to raise its coefficient. Since reinforced plastics are now 
listed among the strongest known materials of construction, this seems like a 
good possibility. When glass fibers are used as the reinforcement, the 
thermal coefficient of expansion is much less th an that of the basic plastic, 
while the stress relaxation (or creep) is considerably decreased. Moreover, 
the strength of the plastic is considerably enhanced. 
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64.2 Radiation Resistance 


The question of radiation resistance of seal materials may arise since 
the vehicles carrying these seals may pass through belts of high intensity 
radiation such as the Van Allen belt. However, the seals would be expected 
to be shielded from direct radiation by other materials, the exact thickness 
and composition of which would be determined by the design. All we can do 
now is outline roughly the behavior which might be expected with polymers. 

The situation is further complicated by the different intensities of radiation 
and different particles encountered in various parts of the radiation belts. 
According to a recent report by R 0 S, Rochlin (Ref. 3), the Van Allen radiation 
consists of protons and electrons trapped by the earth’s magnetic field. The 
protons have energies ranging as high as 700 M.e.v, or more and are capable 
of penetrating up to several inches of lead. The energy distribution and 
number of electrons varies with both positions in the belt and time, but they will 
be almost completely absorbed by thin layers of shielding. For example, a 
detector in the Pioneer IV space probe which was shielded by 1 gram/cm. 2 
(equivalent to a 0.050 inch thickness of stainless steel) indicated a total 
dosage of 10 r./hour in the heart of the inner Van Allen zone. In the heart 
of the outer zones, the dose ranged from 10 r./hour to 100 r. /hour equivalent. 

It was further reduced by a factor of 50 by 4.6 grams/cm, ^ additional shield- 
ing. These doses correspond to 8.3 Rad. /hour maximum in the inner zone and 
83 Rad. /hour maximum in the outer zones. Examination of a report by C.G. 

Collins and V.P. Calkins (Ref. 3) on radiation damage in polymers shows that 
there is no significant damage to most polymers below a dosage of 10^ Rad. ; 
therefore, the effect of Van Allen radiation on polymeric seals can be safely 
ignored in the case of vehicles passing through the belt only once. In the 
case of satellites which may spend significant amounts of time in high 
radiation intensity zones it might be necessary to design seals with some 
radiation resistance. Collins and Calkins list the dosage for threshold 
damage and for 25% damage for a number of polymers, part of which is repro- 
duced on the next page. Since average radiation doses encountered may often 
be less than 10 Rads . /hour and since seals may readily be protected by more 
than the equivalent of 0.05 inches of stainless steel, it is unlikely that 
many cases will arise where radiation resistance will be a prime consideration 
in design. 
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FIGURE 64. 1 


Radiation Damage in Plastics and Rubbers (Ref. 4) 


Material 


Threshold Damage 
Rads . 


25% Damage 
Rads . 


Natural Rubber 

Butyl Rubber 

Neoprene 

Thiokol ST 

Silicone Rubber 

Cellulose Nitrate 

Ethyl Cellulose 

Furan (Dura Ion) 

Asbestos & Carbon Filled 

Methyl methacrylate 
(Plexiglas) 

Trif luorochloroethylene 
(Fluorothene) 

Phenol Formaldehyde 
(Unfilled) 

Polyamide (Nylon FM- 1) 

Polyethylene Terephthalate 
(Mylar) 

Polystyrene 

Polytetraf luoroe thylene 
(Teflon) 

Polyvinyl Chloride 
(Geon 2046) 

Vinylidene Chloride 
(Saran) 


2 x 

10 

6 

2. 5 

X 

10 

2 x 

10 

6 

4 

X 

10 

2 x 

10 

6 

5.5 

X 

10 

5 x 

10 

5 

1.5 

X 

10 

9 x 

10 

5 

5.5 

X 

10 

6.3 

X 

10 5 

5.7 

X 

10 

1.4 

X 

10 6 

5.4 

X 

10 

3.3 

X 

10 8 

3.3 

X 

10 

8.2 

X 

10 5 

1 . 1 

X 

10 

1.3 

X 

10 6 

2.0 

X 

10 

2.7 

X 

10 6 

1. 1 

X 

10 

8.6 

X 

10 5 

4.7 

X 

10 

3.0 

X 

io 7 

1.2 

X 

10 

8.0 

X 

io 8 

4 

X 

10 

1.7 

X 

10 4 

3. 7 

X 

10 

1.9 

X 

io 7 

1 . 1 

X 

10 

4. 1 

X 

io 6 

4.5 

X 

10 


7 

6 

6 

6 

6 

6 

6 

9 

7 

7 

7 
6 

8 
9 

4 

8 

7 


Note: Dosage are defined as follows: 

Threshold dosage: The dosage at which physical property changes first 

become apparent. This represents a point at which the functional 
ability of the material is unchanged. 

25% damage dosage: The dosage at which at least one of the physical 

properties of the material has changed by 25% from its initial value. 
This was arbitrarily chosen as the point beyond which the usefulness 
of the solid materials would be questionable. 

1 1 1 

64-6 



64. 3 Vacuum Effect 

The questions of the effect of a high- vacuum environment on the performance 
of polymeric seals may also be raised. In general, a high vacuum may be 
expected to cause volatilization of low-molecular- weight materials if any 
are present. However, volatilization rate is proportional to vapor pressure, 
and vapor pressure of polymeric materials goes down very fast as molecular 
goes up. Another way of looking at this is to observe that boiling 
point goes up fast with molecular weight. For example, in the linear 
dimethyl siloxane polymer series, (methyl silicones), the boiling point of 
the dimer with a molecular weight of 162 is 100°C at 1 atmosphere. As the 
series progresses to longer and longer chains, the boiling point rises very 
fast. At molecular weight 278, it is 153°C, at molecular weight 394, it 
is 194°^ and at a molecular weight of 742, the boiling point is 199 C at a 
pressure of 0.02 atmospheres. Since rubbers have an infinite molecular 
weight because they are crosslinked, and even uncrosslinked straight chain 
polymers have molecular weights in the hundreds of thousands and in the 
millions, they have substantially zero vapor pressure; therefore, they can 
be regarded as non-volatile. Some plastics do, however, contain plasticizers 
which are lower- molecular- weight materials added to increase flexibility. It 
would be possible to distill out some of these molecules in long periods under 
high vacuum, particularly at high temperature. The best way to avoid this 
possibility and consequent property changes is to choose materials without 
low-molecular-weight plasticizers for use where long exposure under high 
vacuum is necessary. It is usually necessary to confer with the supplier to 
make this choice, since it is impossible to tell from trade names which 
formulations contain low-molecular-weight plasticizers. 



64. 4 Oxidation Degradation 


One effect of high-vacuum environment will be a benefit due to the 
absence of oxygen. The most rapid deterioration of many plastic materials 
occurs by a mechanism of oxidation. This effect has been shown to be due to 
the oxygen in the air which chemically attacks the rubber molecules, cutting 
them like a pair of scissors. This has been observed by a host of workers, 
Refs. 8 to 10, who have found that the experimental approach of stretching 
the rubber in an air oven and measuring the decrease in force required to 
hold the same extension is a convenient way of measuring the rate of attack 
by oxygen. Sometimes the oxygen causes crosslinks to form between chains 
faster than the degradative mechanism can cut chains. In this case the poly- 
mer may stiffen with time in oxygen at high temperature rather than soften, 
as is the case when, the chain cutting or scission reaction predominates. In 
either case, the end product is a deterioration in properties from the stand- 
point of the seal designer. If scission predominates, the seal will gradually 
loosen with time as the stress relaxes. If crosslinking predominates, the 
rubber will become too hard and lose its flexibility, so that it may even^ 
tually crack. These processes have been so well documented for some materials 
that the stress-relaxation curves for different temperatures can be shown to 
obey the Arrhenius rate law. Figure 64,2 shows how the stress relaxation due 
to oxidation of natural rubber follows the same curve for different tempera- 
tures but displaced along the logarithmic time scale. Such a displacement 
indicates a process following a simple rate law and having a constant acti- 
vation energy. 
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FIGURE 64. 3 


STRESS RELAXATION OF POLYSULFIDE RUBBER AT VARIOUS ELONGATIONS AT 60°C (Ref. 6) 

(Stress -jr initial stress) 



• 10% Elongation 
O 30% Elongation 
© 50% Elongation 
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64.. 5 Thermal Degradation 


Naturally, none of these effects will be observed in the absence of 
oxygen, but there are other types of degradation which can occur in polymers 
even without oxygen. If the temperature is high enough, polymer chains 
may begin to split up into smaller fragments, which may volatilize if 
small enough. Sometimes this process occurs in a kind of polymerization 
in reverse or chain unzipping. When this happens, polymer molecules 
lose one monomer unit at a time like beads coming off a string. The 
monomers are quite volatile at the temperature where this occurs, so they 
boil off immediately. In polystyrene, for example, this process begins 
at 300°C. 

Another type of degradation at elevated temperature is bond exchange. 
This results in stress relaxation just as bond scission by oxidation does, 
but the chemistry is a little different and it does not depend on the 
presence of oxygen. In polysulfide rubber, for example, bond exchange 
occurs between sulfur atoms on adjacent Q chains , resulting in rapid stress 
relaxation at temperatures as low as 60 C, as illustrated in Fig- 64.3. 


A convenient way of ranking the thermal stability of polymers is the 
use of thermogravimetric analysis. This test consists in heating the 
polymer sample at a constant rate while accurately measuring its weight 
change. This is most often done in nitrogen to avoid the effects of 
oxidation. The weight curve will show a sharp drop where decomposition 
starts and the temperature at which this occurs can be regarded as a measure 
of the thermal stability of the polymer. 

FIGURE 64.4 


Decomposition Temperatures of 
Some Common Polymers as Determined by 
Differential Thermal Analysis 


Material Temperature , C 


Polystyrene 

395 

Maleic hardened epoxy 

405 

Polymethylmethacrylate 

345 

6-6 Nylon 

419 

Polytetraflu or ethylene (Teflon) 

555 

Polytr if luorochloroethylene (Kel -F) 

410 

Viton A 

460 

Silicone resin 

505 
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64. 6 Differential Expansion 

As mentioned briefly earlier, the differential coefficient of expansion 
between polymers and metals will result in the loosening of a seal compressed 
by metal elements as the temperature is lowered below the glass transition. 
Furthermore, it will make almost no difference how much the seal is compressed 
when it is in the rubbery state; as the polymer cools through the glass 
transition, chain coiling and uncoiling is frozen and the only motions left 
to the polymer molecules are small displacement of the chain segments with 
respect to other segments on the same molecule and segments on adjacent 
molecules. As the polymer is cooled, it would normally shrink, but if it 
is restrained from doing this by a metal framework, the chain segments will 
be effectively displaced from their equilibrium positions by an amount equal 
to the differential coefficient of expansion between the metal and the 
polymer. Therefore, the force exerted by the polymer molecule on its 
restraining elements can be calculated roughly by using the strain due to 
the differential coefficient of expansion and multiplying by the average 
Young's modulus of elasticity over the temperature interval considered. 

An example of this behavior is shown in Fig. 64.5. These data were obtained 
for polymethylmethacrylate (Plexiglas) which exhibits typical rubber— like 
behavior at temperatures above 105°C and has a glass transition between 
80°C and 105°C. The experiments illustrated show the change in force 
required to hold fixed the elongation of a sample of this material while 
lowering the temperature. Two different elongations were used;^in each 
case the elongation was applied by stretching the sample at 150 C and 
then the oven temperature was allowed to drog slowly to room temperature. 

As soon as the temperature dropped below 100 C, which is in the glass 
transition, the force (plotted as stress) began to rise and continued 
steadily up until the oven reached the lowest temperature. The oven was 
then warmed and the stress fell almost on the same curve with some hysteresis 
at the glass transition. The important observation is that with two widely 
different elongations imposed at 150°C, the identical curve was traced below 
the glass transition. If one assumes that all the strain below the glass 
transition is due to linear coefficient of expansion of the polymer 
(9 x 10" 5 /°C) and, using an average Young's Modulus of 300,000 P.S.I^ in 
the glassy region, one calculates a stress of 1350 P.S.I. or 93 x 10 
dynes /cm. which is almost exactly what is observed. Thus, the polymer 
sample "stretched itself" as it cooled. Furthermore, if it had been in a 
compressed state upon cooling through the glass transition, it would have 
made little difference; it still would have contracted and either stretched 
itself or pulled away from its restraining elements unless it was fastened 
to them by an adhesive. 
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64. 7 Glass Transition 


The behavior exhibited by polymethylmethacrylate in the above illustration 
is a perfectly general one, which will be exhibited by any polymer upon 
cooling through the glass transition; and it illustrates the importance of 
knowing the glass transition temperature, the coefficient of expansion of 
the polymer, and its modulus of elasticity. A further observation in the 
above example is how closely the stress levels on the warming curve agree 
with the stress levels on the cooling curve; this illustrates the relatively 
low degree of stress relaxation which has occurred in the glassy state due to 
the progressive freezing-in of modes of relaxation as the temperature is 
lowered. 

Since the glass transition is so important in the behavior of polymers, 
a short table is reproduced below to show where this transition falls for 
some common polymers. It will be noted that all the polymers listed encounter 
this transition upon cooling to -100°C. 


FIGURE 64. 6 

Polymer Glass Transition Temperatures 


Material Glass Transition, °C 


Poly isobutylene 

-74 

Natural Rubber (Hevea) 

-73 

Polyvinyl acetate 

+32 

Polyurethane 

-35 

Polystyrene 

100 

Polymethylmethacrylate (Plexiglas) 

105 

Polyvinyl chloride 

74 

Butadiene - styrene rubber 

-33 

Silicone rubber 

-80 

Polytetraf luoroethylene (Teflon) 

+ 25 


The way in which the modulus of elasticity increases sharply as temperature 
decreases below the glass transition is illustrated in the curves shown in 
Fig. 64.7. These depict the shear modulus of several rubbers as a function 
of temperature. It will be noted that an increase in modulus ranging between 
a factor of 100 and 1000 occurs in a fairly narrow temperature interval. 
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FIGURE 64.7 

YOUNG'S MODULUS VS. TEMPERATURE FOR VARIOUS VULCANIZED RUBBERS (Ref. 6) 

O natural rubber 
y polybutadiene 
y butadiene-styrene 85/15 
□ butadiene -styrene 75/25 
A butadiene-styrene 60/40 



64. 8 Glassy State Stress Relaxation 


The properties of polymers in the glassy state below the glass 
transition are an important consideration where polymers are to be used as 
seals at low temperature. Some polymers have been studied directly in this 
temperature region by use of stress relaxation measurements (Refs, 11, 12). 
These measurements show that stress relaxation does indeed occur in glassy 
polymers, but at a very slow rate compared to relaxation in the transition 
region. The relaxation is, of course, a function of temperature and depends 
on free volume and often on water content and plasticizer content. The free 
volume in a polymer below its glass transition depends on thermal history. 
Samples cooled rapidly have greater free volume than those cooled more 
slowly; therefore, rapidly cooled samples exhibit greater relaxation of stress 
at low temperature than those which have more time to accommodate to the 
temperature change during the cooling process. 

Because relaxation is a slow process in the glassy state, and due to 
experimental difficulties in measurements at low temperatures over long 
periods of time, few relaxation measurements have been made directly on poly- 
mers below the glass transition where the glass transition occurs at low 
temperature. Instead, the technique of dynamic measurements has been widely 
used to study the response of polymers at small deflections below their glass 
transition temperatures. This technique consists in subjecting the specimen 
to a cyclic load and measuring the amplitude and phase shift of the resulting 
deflection. The ratio of stress amplitude to strain amplitude is the complex 
modulus of elasticity, and the phase angle between the stress amplitude and 
strain amplitude is the loss angle 5. For a material with zero loss, the 
stress amplitude and strain amplitude would be in phase, and the angle 5 would 
be zero. The tangent of the angle 5 indicates the relative magnitudes of 
viscous and elastic effects in the material, and hence the relative degree of 
stress relaxation to be expected. For a simple Maxwell (spring - dashpot) 
viscoelastic model of the material, the relaxation time is the time for 63% 
decay in stress under a steady-state deformation and is 10/n times the period 
of oscillation for which tan 5 = 0.05. 


By making measurements at different frequencies and temperatures, it is 
sometimes possible to estimate long-term behavior at low temperatures from 
higher- frequency behavior at higher temperatures. In order for this to be 
done, it is necessary that all parts of the frequency- response curve have the 
same temperature dependency. This is a fair assumption if all the measure- 
ments are made far enough below the glass transition. At the present time, 
however, only limited frequency and temperature response data are available 
for most polymers at low temperatures, so only qualitative remarks concern- 
ing their viscoelastic behavior can be made. That is, it can certainly be 
said that polymers which exhibit low dynamic losses below their glass transi- 
tions will exhibit low stress relaxation. Some dynamic measurements on 
Plexiglas and polyethylene are shown in Figures 64.8 and 64. 9 , respect ively. 
Examination of some of these dynamic- losses- versus- temperature plots which 
contain data at several frequencies shows that the loss curves shift to lower 
temperature in a fairly uniform manner as the frequency is reduced. For 


64-16 


121 



example, in Figure 64,8, the point where tan 5 is 0.05 shifts from 35°c to 
-20°C as the frequency is reduced from 2000 c.p.s. to 0.4 c.p.s. If one 
assumes that the shift is uniform on a logarithm-of- frequency scale, the 
temperature at which tan & is 0.05 at a frequency of 8 x 10"^ c.p.s. would be 
-75°C. A frequency of 8 x 10“ ^ c a p,s. corresponds to a relaxation time of 
40,000 seconds, or about half a day. This means that there would be sub- 
stantial decay in stress observed in Plexiglas when strained at -75 C and 
held at constant strain during a period of half a day. However, such a stress 
relaxation on this time scale would fall off fast at lower temperatures, 
since the loss is decreasing with temperature at the frequency for which tan 
5 = 0,05. These are only qualitative interpretations of the data, which is 
all that is possible where the frequency- loss curves do not slide perfectly 
along the temperature scale, as evidenced by the fact that the loss peak 
height is much greater at 2000 c.p.s, than at 0.4 c.p.s. The fact that the 
loss peak height decreases with reduction in temperature would suggest that 
the relaxation would be even lower at -75 C than the above simple extra- 
polation might indicate. 
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Temperature 

FIGURE 64.8 

MECHANICAL DAMPING OF POLYMETHYLMETHACRYLATE 
(PLEXIGLAS) AS A FUNCTION OF TEMPERATURE AND FREQUENCY (Ref.l) 
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FIGURE 64.9 


MODULUS AND LOSS ANGLE TANGENT , TAN £ , OF UNORIENTED CRYSTALLINE 
POLYETHYLENE TEREFHTHALATE AT VARIOUS FREQUENCIES (Ref. 1) 


0 10,000 c.p.s. 
<8 5.9 c.p.s. 

X 0.72 c.p.s. 

O 0.082 c.p.s. 
^ 0.0094 c.p.s. 
□ 0.001 c.p.s. 


12 ? 
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64,9 Ultimate Strength and Elongation at Low Temperature 


When a limited amount of dynamic measurements are available at too few 
frequencies, it is difficult to carry out the necessary correlations between 
dynamic measurements and stress relaxation. More measurements of dynamic 
properties at low temperature are needed to predict low- temperature stress 
relaxation, but an even greater need is for ultimate- strength measurements on 
polymers at low temperatures, since these materials can be expected to be 
brittle below the glass transition temperature; the lower the temperature 
below the glass transition, the greater the likelihood of brittleness will.be 
Yet, only limited data are available (Ref. 13) to show either how much stress 
various polymers will endure at very low temperatures or how much elongation 
they will endure before breaking. These measurements are probably of greater 
importance than stress relaxation, creep, or dynamic measurements since fail- 
ure by cracking at these temperatures seems to be a much greater possibility 
than failure by distortion. It must be remembered that no known rubbers act 
like rubbers at liquid-oxygen and liquid-nitrogen temperatures. All are more 
or less brittle, glass-like materials, just as are all plastics. Which 
materials are least subject to low elongation failure would seem to be the 
key question, and this can only be answered by experiment. 
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64, 10 Conclusions 


Exposure to the intensity and kind of radiation exhibited in the Van 
Allen belts will not have a serious effect on polymer properties unless 
exposure times of years are required or unless the polymeric materials 
are exposed without shielding. 

The effect of high vacuum on polymer properties will be negligible 
unless the formulation contains more than small amounts of low molecular 
weight plasticizers or solvents. 

High temperatures cannot be tolerated by polymers very long. Many 
rubbers oxidize and either soften or embrittle at temperatures of 150°C 
(302°F) and lower. Others undergo chemical changes without oxygen. 
Thermoplastics such as polystyrene and polymethylmethacrylate revert to 
monomer at temperatures over 300°C (572°F) and their properties deteriorate 
slowly at much lower temperature. Polysulfides deteriorate at 60 C (140 F). 
Higher temperature polymers are in development but 250-300°C looks like the 
upper operation limit at present for these materials. 

Low temperature causes all polymeric materials to stiffen because 
they pass through glass transitions before reaching -100 C (-148°F). In 
their glassy state below this transition, all polymers are very stiff and 
more or less brittle although very little information is available on 
strengths or failure elongations at liquid oxygen and liquid nitrogen 
temperatures. Dynamic measurements havq been made at such low temperatures; 
these reveal high moduli and very little flow. 

All polymers have much greater thermal expansion coefficients than 
metals and, therefore, present problems if used in contact with metals at 
low temperature where the polymers have little flexibility or flow and, 
hence, cannot accommodate as well. The use of fillers to decrease expansion 
coefficients of polymers is a possibility, but most fillers increase stiffness 
and decrease the failure elongation in the glassy state. The use of fiber- 
reinforced, polymeric supporting elements in contact with polymeric gaskets 
may help relieve the problems due to coef ficient-of-expansion differences. 
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